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Abstract 

We report on determinations of \V ub \ resulting from studies of the branching fraction and q 2 
distributions in exclusive semileptonic B decays that proceed via the b — > u transition. Our 
data set consists of the 9.7 x 10 6 BB meson pairs collected at the T(45) resonance with the 
CLEO II detector. We measure B(B° -► tt~£ + u) = (1.33 ± 0.18 ± 0.11 ± 0.01 ± 0.07) x 10~ 4 
and B(B° -> p~£ + v) = (2.17 ± 0.34 +g;f] ± 0.41 ± 0.01) x 10~ 4 , where the errors are statistical, 
experimental systematic, systematic due to residual form-factor uncertainties in the signal, and 
systematic due to residual form-factor uncertainties in the cross-feed modes, respectively. We 
also find B(B + -► -qt+v) = (0.84 ± 0.31 ± 0.16 ± 0.09) x 10~ 4 , consistent with what is expected 
from the B — > ir£v mode and quark model symmetries. We extract \V u b\ using Light-Cone Sum 
Rules (LCSR) for < q 2 < 16 GeV 2 and Lattice QCD (LQCD) for 16 GeV 2 < q 2 < q 2 nax . 
Combining both intervals yields \V ub \ = (3.24 ± 0.22 ± 0.13 IqIq ± °- 09 ) x 10 ~ 3 for and 
\V uh \ = (3.00 ±0.21 toll -0.38 ±0- 28 ) x for p£u, where the errors are statistical, experimental 
systematic, theoretical, and signal form-factor shape, respectively. Our combined value from both 
decay modes is \V ub \ = (3.17 ± 0.17 t° f 7 ±8;|| ± 0.03) x lO" 3 . 
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I. INTRODUCTION 

The element V u b remains one of the most poorly constrained parameters of the Cabbibo- 
Kobayashi-Maskawa (CKM) matrix [1]. Its magnitude, \V u b\, plays a central role in con- 
straints based on the unitarity of the CKM matrix and inputs from both CP-conserving 
processes in the B meson decay and CP-violating processes in the neutral kaon and B 
systems. The value of \V u b\ and, in particular, the accuracy to which we have measured this 
important parameter, have been the subjects of considerable debate over the past decade 
[2]. An accurate determination of \V u b\ with well- understood uncertainties remains one of 
the fundamental priorities for heavy flavor physics. 

A number of \V u b\ measurement approaches have been attempted, and are reviewed in 
reference [2]. Inclusive techniques are hampered by a mismatch in kinematic regions where 
the large experimental backgrounds from b — > civ can be suppressed versus regions in which 
the theoretical uncertainties can be reliably determined. For exclusive reconstruction of 
particular final states, the primary challenge is calculation of the form factors for those 
channels. The first measurements of exclusive charmless semileptonic branching fractions 
[3], including evaluation of \V u b\, were performed by the CLEO experiment at the Cornell 
Electron Storage Ring (CESR) using the modes B° -> tt-£ + v, B + -> tt°£ + v, B° -> p~£ + v, 
B + — > p°£ + u, B + — > uj£ + u, and charge-conjugate decays, where £ = e or p. A second 
measurement of the p£v modes by CLEO [4], using similar techniques but a much different 
signal to background optimization, provided consistent, essentially independent, results with 
a similar total uncertainty. The combined analyses yielded \V ub \ = (3.25±0.14tg;29 ±0.55) x 
10~ 3 , where the errors are statistical, experimental systematic and estimated theoretical 
uncertainties, respectively. The ir and p modes contribute about equally to this result. 

This paper presents an update of the original exclusive B — > X u £u analysis [3], and is 
based on a total data sample of 9.7 x 10 6 BB pairs collected on the T(AS) resonance. The 
results presented here supersede those of reference [3]. In addition to using a larger data 
set, the analysis has been modified to minimize uncertainties arising from the momentum- 
transfer (q 2 ) dependence of the form factors. Most notably, the lower bounds on the charged- 
lepton momentum for both the pseudoscalar and the vector modes have been lowered, and 
the branching fractions are determined independently in three q 2 regions. For the p modes, 
the branching fractions as a function of q 2 were first determined by the second CLEO 
p£v analysis [4]. The present analysis has a significantly broader accepted range for the 
charged lepton momentum, which allows for better discrimination among models. A detailed 
description of this analysis can be found in reference [5]. 

II. EXCLUSIVE CHARMLESS SEMILEPTONIC DECAYS 

The semileptonic transition of a B meson (a pseudoscalar) to a final state with a single 
pseudoscalar meson P can, in the limit of a massless charged lepton, be described by a single 
form factor fi(q 2 ): 



where y = q 2 /Mg, M B is the mass of the B meson, Gf is the Fermi constant, kp is the 
meson momentum, and 9we is the angle between the charged lepton direction in the virtual 
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W {£+v) rest frame and the direction of the virtual W in the B rest frame. For a transition to 
a final state with a single vector meson V, three form factors (Ai, A 2 , and V) are necessary: 



- \ Vub \ mo_3 X l Z J 



dy d cos W e 128n 3 

[1 — cosO W e) 2 } n hl + (1 + cosfl^i) 2 1 o' + s i n2 ^^l-^t ' 2 



\HJ 2 , „ ,J#J 2 



where ky is the meson momentum and the three helicity amplitudes are given by 

1 
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M B + m v 
1 M B 



A 1 (q 2 )T2M B k v V(q 2 ) , and (3) 

(4) 



^y2m v {M B + my) [\ M 2 B 



The structure of these differential decay rates immediately allows us to draw some general 
conclusions regarding the properties of the semileptonic decays that we reconstruct in this 
analysis. For the p(uj)iv transitions, the left-handed, V — A, nature of the charged current 
at the quark level manifests itself at the hadronic level as \H-\ > \H + \. The H- contribution 
is also expected to dominate the Hq contribution, leading to a forward-peaked distribution 
for cos 9we- For n(r))£v, there is a sin 2 Owe dependence, independent of the form factor. 
The pseudoscalar modes also contain an extra factor of the meson momentum squared, 
which suppresses the rate near g 2 ^ {hp = 0). Taken together, these two effects give the 
pseudoscalar modes a softer charged lepton momentum spectrum than the vector modes. 

Calculation of the form factors has become a considerable theoretical industry, with a 
variety of techniques now being employed. Form factors based on lattice QCD (LQCD) cal- 
culations [6-18] and on light-cone sum rules (LCSR) [19-27] currently have uncertainties in 
the 15% to 20% range. A variety of quark-model calculations exist [28-42]. Finally, a num- 
ber of other approaches [43-48], such as dispersive bounds and experimentally-constrained 
models based on heavy quark symmetry, all seek to improve the range of q 2 over which 
the form factors can be estimated without introduction of significant model dependence. 
Figure 1 illustrates the broad variation in shape that arises in the literature. Unfortunately, 
all the form-factor calculations currently have contributions to the uncertainty that are un- 
controlled. The light-cone sum rules calculations assume quark-hadron duality, offering a 
"canonical" contribution to the uncertainty of 10%, but with no known means of rigorously 
estimating that uncertainty. The LQCD calculations to date remain in the "quenched" ap- 
proximation (no light quark loops in the propagators), which limits the ultimate precision 
to the 15% to 20% range. With the quark-model calculations it is difficult to quantify the 
uncertainty of a particular calculation by their very nature. These uncertainties in the form 
factors translate directly into the same fractional uncertainty on \V u b\. 

In the Tciu modes, with only a single form factor in the massless lepton approximation, 
we expect that the rates extracted in the q 2 intervals that we have chosen will be largely 
independent of the form-factor shapes. In the vector modes, however, the three form factors 
interfere and differences in this interference among models, particularly at lower q 2 values, 
can lead to a residual model dependence. To investigate this effect, we will analyze the 
vector modes with three separate charged lepton momentum requirements. 
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FIG. 1: Predictions for dT(B — » Tttv)jdq 2 (left) and for dT{B — > plu)/dq 2 (right) for a variety of 
calculations, illustrating the range of variation of the predicted g 2 -dependence. See Section VI for 
further discussion of the calculations. 



III. EVENT RECONSTRUCTION AND SELECTION 

The CLEO detector [49, 50] contains three concentric tracking devices within a 1.5 T 
superconducting solenoid that detect charged particles over 95% (93%) of the solid angle 
for the first third (last two thirds) of the data. For the last two thirds of the data, a silicon 
vertex detector replaced a straw-tube wire chamber. The momentum resolution at 2 GeV/c 
is 0.6%. A CsI(Tl) electromagnetic calorimeter, also inside the solenoid, covers 98% of 4tt. 
A typical tt° mass resolution is 6 MeV. Charged tracks are assigned the most probable mass 
based on specific ionization, time of flight, and the relative rates as a function of momentum 
for proton, K + , and 7r + production in B decay. 

The undetected neutrino complicates analysis of semileptonic decays. Because of the 
good hermeticity of the CLEO detector, we can reconstruct the neutrino via the missing 
energy (E m i SS = 2i?beam — J2Ei) and missing momentum (P m i SS = —J2f>i) m each event. 
In the process e + e~ — > T(45) — ► BB, the total energy of the beams is imparted to the 
BB system; at CESR, that system is at, or nearly at, rest. (A small crossing angle has 
been in use at CESR for most of the running.) The missing mass, M^ iss = E^ iss — |-P m i SS | 2 , 
must be consistent, within resolution, with a massless neutrino. Specifically, we require 
-0.5 < M^ iss /2E miss < 0.3 GeV for events with a total charge AQ = 0, and |M^ iss |/2E miss < 
0.3 GeV for events with | AQ\ = 1. 

Signal Monte Carlo (MC) events show a |P m iss| resolution of 85 MeV/c. The resolution on 
-E m i SS is about three times larger than the momentum resolution [51]. Significant effort has 
been devoted to minimizing multiple counting of charged particles in the track reconstruction 
(e.g., particles that curl multiple times within the tracking volume), and to suppressing 
clusters in the calorimeter from charged hadrons that have interacted. 

With an estimate of the neutrino four-momentum in hand, we can employ full recon- 
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struction of our signal modes. Because the resolution on -E m i SS is so much larger than 
that for |P m i ss |, we use (E v ,p v ) = (]P m i ss |, -P m is S ) for full reconstruction. The neutrino 
combined with the signal charged lepton (£) and meson (m) should satisfy, within reso- 
lution, the constraints on energy, AE = [E v + E e + E m ) — E hcam m 0, and on momentum, 
M m i v = [E 2 c&m — \ a Vv + Pi + Pm| 2 ] 5 ~ M B , where a is chosen to force AE = 0. The neutrino 
momentum resolution dominates the AE resolution, so the momentum scaling corrects for 
the mismeasurement of the magnitude of the neutrino momentum in the M m £ U calculation. 
Uncertainty in the neutrino direction then remains as the dominant source of smearing in 
this mass calculation. 

We reconstruct q 2 = M^* = (p u + Pi) 2 f° r eac h decay from the reconstructed charged 
lepton four-momentum and the missing momentum. In addition to using the scaled recon- 
structed momentum ap described above, the direction of the missing momentum is changed 
through the smallest angle consistent with forcing M mtu = M B . This procedure results in 
a q 2 resolution of 0.3 GeV 2 , independent of q 2 . The irlv and the plv modes are analyzed 
separately in the intervals q 2 < 8 GeV 2 , 8 < q 2 < 16 GeV 2 , and q 2 > 16 GeV 2 . For the uj£u 
and r\iv modes, for which we have low statistics, we sum over all q 2 . 

Information from specific ionization is combined with calorimetric and tracking measure- 
ments to identify electrons with p t > 600 MeV/c over 90% of the solid angle. Particles 
registering hits in counters deeper than 5 interaction lengths over the polar angle range 
| cos 6*| < 0.85 are considered muons. Those with hits beyond 3 interaction lengths over 
| cos#| < 0.71 are used in a multiple-lepton veto, described below. Candidate leptons must 
have pt> 1.0 GeV/c for the n and rj (pseudoscalar) modes, and pt > 1.5 GeV/c for the p 
and uj (vector) modes, which can couple to the W helicities ±1 and hence have a harder 
spectrum. This momentum requirement for the vector modes defines the nominal analysis. 
We also analyze the vector modes with the lepton momentum requirements pu > 1.75 GeV/c 
and pi > 2.0 GeV/c. The identification efficiency above 1.5 GeV/c averages over 90%; the 
probability that a hadron is misidentified as an electron (muon), a fake lepton, is about 0.1% 
(1%). 

The 5-interaction-length requirement for muons causes the muon acceptance to fall 
rapidly below 1.4 GeV/c. As a result, only electrons contribute at the low end of the 
momentum range we accept for niu, and electrons dominate the measurement in the lowest 
q 2 interval. 

A 7T° candidate must have a 77 mass within 2 standard deviations of the n° mass. We 
reconstruct the u via its 7r + 7r~7r° decay, reducing combinatoric background by rejecting 
combinations away from the center of the ui Dalitz plot. We reconstruct r\ in both the 77 
and the 7r + 7r~7r° decay modes. For the 77, we require the reconstructed mass to be within 
2 standard deviations of the 77 mass (within about 26 MeV). For the 7r + 7r~7r°, we require 
IrriK+K-nO — m v \ < 10 MeV (about 1.7 times the resolution). We impose a kinematic mass 
constraint on the momentum of all 7r° or rj candidates in the 77 final state. 

Backgrounds arise from the e + e _ — > qq and e + e~ — > t + t~ continuum, fake leptons, 
b — > civ, and B — > X u lv modes other than the signal modes. Backgrounds from continuum 
processes are suppressed by use of two event-shape variables. The selection criteria were 
optimized using background and signal Monte Carlo samples, rather than data, to avoid 
potential bias. The first variable is the angle (cos thrust) between the thrust axis evaluated 
for the candidate signal-mode particles (not including the neutrino) and that for the rest 
of the event. (The thrust axes are signed by picking the hemisphere containing the most 
energy.) For BB events at CESR, the distribution in this variable is flat because the £>'s 
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are nearly at rest and thus their decay orientations are independent. For continuum events 
the distribution is strongly forward and backward peaked. The ratio R2 of the second to the 
zeroth Fox- Wolfram moment [52], which distinguishes spherical from jetty topologies, is also 
utilized. The continuum background tends to have a small reconstructed q 2 . We therefore 
tune the continuum cut employed in the i? 2 -cos thrust plane separately in each q 2 interval, 
and separately for the n and p modes. Signal events with low q 2 appear rather jetty, so 
a cut using R 2 , when data is binned over a broad q 2 range, would introduce an efficiency 
bias. So for the uj and rj modes, for which all q 2 regions are combined, only a cos thrust cut 
is applied, reducing uncertainties from the g 2 -dependence of the form factors. Our criteria 
suppress the continuum background by over a factor of 10 and are about 80% efficient. 

The \pi\ cuts greatly reduce background from b — > c — > siv and bias mildly against 
b — > civ. For the vector modes, we further require cos 6^ > 0, since the signal rate is 
largely suppressed by V — A outside this region, while the background is roughly flat in the 
region excluded, and falls off in the region accepted. 

Backgrounds, particularly b — > civ, can smear into the signal region in AE and M m p v 
when -P m i SS misrepresents p u . Such backgrounds are highly suppressed by rejecting events 
with multiple charged leptons or a total event charge \AQ\ > 1, both of which indicate 
missing particles. Requiring M^ iss to be consistent with zero also provides powerful back- 
ground suppression. Still, Monte Carlo studies show that the dominant remaining b — > civ 
events contain either a Kl meson or a second neutrino (from c — > siv, with the lepton not 
identified) that is roughly collinear with the primary neutrino. 

Our selection criteria studies, based on statistical considerations, indicated that keeping 
the \AQ\ = 1 sample as well as the AQ = was favorable in spite of the poorer signal-to- 
background ratio. Further systematic considerations indicated that the use of the \AQ\ = 1 
sample remained advantageous for the pseudoscalar modes. For the vector, in particular the 
p modes, however, the overall poorer signal-to-background ratio made the \AQ\ = 1 sample 
overly sensitive to systematic effects in both the modelling of the B — > X u iv backgrounds 
and the simulation of the detector. Therefore for the vector modes we require AQ = 0. 

IV. EXTRACTION OF BRANCHING FRACTIONS 
A. Method and binning 

To extract the branching fraction information, we performed a binned maximum likeli- 
hood fit that was extended to include the finite statistics of the Monte Carlo, off-resonance, 
and fake-lepton samples following the method of Barlow and Beeston [59]. The data 
in each mode were coarsely binned over the two dimensional region 5.175 < M m £ U < 
5.2875 GeV, \AE\ < 0.75 GeV. We further binned the data in the reconstructed 2tt and 
37r masses in the p and uj modes. The |AQ| = 1 samples were binned separately from 
AQ = samples. Separation of the net charge samples allowed us to take advantage of the 
better signal-to-noise ratio of the \AQ\ = sample while reducing our dependence on our 
knowledge of the absolute tracking efficiency Finally, we binned the data in q 2 for the two 
ixiv and the two piv modes. For the ujIv and the rjiv modes, we combined all q 2 information 
into a single bin. 

Our fitting strategy was designed to minimize dependence of the results on the details of 
the simulation - both from detector and physics standpoints. The choice of binning balanced 
separation of signal and background against reliance on detailed MC shape predictions. To 
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TABLE I: Summary of the number of bins used in each mode for the nominal fit. 



AE, M mlv AQ M 2wM g 2 total 

■K-l+v 7 2 1 3 42 

tt°£ + u 7 2 1 3 42 

p~l+v 7 1 3 3 63 

p°l + i> 7 1 3 3 63 

col+u 7 1 3 1 21 

rj^i+it 7 2 1 1 14 

r^t+v 7 2 1 1 14 



help minimize the model dependence of the branching fraction determinations, we did not use 
information from the lepton momentum spectrum or from cos 6we within the fit. Extraction 
of rates in the separate q 2 intervals further reduces reliance on the form factors. 

The AE bin intervals used in the nominal fit were —0.75 < AE < —0.45 GeV, —0.45 < 
AE < -0.15 GeV, and -0.15 < AE < 0.25 GeV (the AE signal band). The M mlv bin 
intervals were 5.175 < M mlv < 5.2425 GeV and 5.2425 < M mlv < 5.2875 GeV. In the AE 
signal band, this second mass interval is divided into two equal bands. Hence we used a 
total of seven bins in these two variables. In the piv {uilv) modes, we used three equal bins 
over the 2ir (3tc) mass range within ±285 MeV (±30 MeV) of the nominal p (u) mass. The 
three q 2 intervals in the iriv and the piv modes were q 2 < 8 GeV 2 , 8 < q 2 < 16 GeV 2 , 
and q 2 > 16 GeV 2 . The number of bins for each mode in the nominal fit is summarized in 
Table I. The nominal fit had a total of 259 bins. For studies in which the \AQ\ = 1 sample 
is included in the p and uo modes, the fit had an additional 147 bins for a total of 406 bins. 

To examine yields, efficiency, and kinematics in this paper, we use the most sensitive 
bin (the "signal bin") 5.265 < M mlu < 5.2875 GeV and -0.15 < AE < 0.25 GeV, though 
neighboring bins also contribute information to the fit. For comparison, the M m £ U and AE 
resolutions are about 7 MeV and 100 MeV, respectively, dominated by the resolution on \p u \. 
The 2ir (or 3n) mass intervals ±95 MeV and ±10 MeV, centered on the nominal masses, 
are used for figures for p and uo candidates, respectively. 

To simplify the statistical interpretation of the results, we limited the number of multiple 
entries per event. For each individual mode, the candidate with the smallest \AE\ among 
those satisfying M mlv > 5.175 GeV was chosen, independent of q 2 . A given event could 
contribute to multiple modes, although contribution near the signal region in more than one 
mode was rare. In the p and uo modes, each of the mass bins described above was considered 
a separate mode. 



B. Fit components and parameters 

MC simulation provided the distributions in each mode for signal, the b — > c background, 
the cross feed among the modes, and the feed down from higher mass B — > X u lv decays. It 
included a full description of the b — > c and charm decay modes and a GEANT-based [53] 
detector model. The X u iv feed down was evaluated with a simulation of the B — > X u iv 
process based on an inclusive operator product expansion (OPE) calculation [54] of dT(B — > 
X u lv)/dMx u , using parameters determined from the CLEO analysis of the B — > X s ^ photon 
spectrum [55, 56] (also used in the recent CLEO lepton-momentum end-point analysis [57]). 
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The nominal analysis combined this inclusive spectrum with the ISGW II model [31] for 
all mesons through the p(1450). For each exclusive mode, we "subtracted rate" from the 
inclusive calculation with a weight of the form exp [—a(Mx u — Mr) 2 / Aq CD ], where Mr is the 
central mass of the resonance R. At any given M Xu , the rate remaining after this subtraction 
of the exclusive modes is hadronized nonresonantly. Variations of the inclusive parameters 
based on the uncertainties in the B — > X s 7 analysis and variations of the hadronization 
model {e.g., fully nonresonant but with Tmiu removed from the p mass region) are included 
in the systematic uncertainties. The signal modes are excluded from these B — > X u iv 
samples. 

The contributions from events in which hadrons have faked the signal leptons and from 
continuum are evaluated using data. The electron and muon identification fake rates from 
pions, kaons, and protons are measured in data using a variety of tagged samples. The 
analysis is performed on a sample of hadronic events with no identified leptons, treating 
each track in turn as a signal electron and then a signal muon. The contribution in each 
mode is weighted according to the fake rate. 

We determined the residual continuum background using data collected 60 MeV below 
the T(4S) energy. The center-of-mass energy and cross-section differences were taken into 
account as necessary. For each combination of mode, reconstructed q 2 bin, and for each AQ 
value, we determined the rate over the full AE — M m £ U plane by applying all cuts, including 
continuum-suppression cuts, and then scaling according to the relative on-resonance and 
off-resonance luminosities. To smooth the statistical fluctuations within each combination, 
we determined the shape over the AE — M miv plane by the following procedure. First, we 
dropped the continuum-suppression cuts, and obtained the shape over the AE — M miv plane, 
for each combination from data. Then, from continuum qq MC, t + t - MC, and our fake 
lepton samples, we determined the change in shape over the AE — M m £ U plane caused by 
application of the continuum-suppression cuts, i.e., we obtained the ratio of yields, with to 
without cuts, for each AE — M m i u bin, for each combination. Applying the ratios so obtained 
to the off-resonance data without continuum- suppression cuts, we obtained the shape of the 
background over the AE — M m i v plane, for each combination. 

For each signal mode, we generated a sample of signal Monte Carlo that is flat in phase 
space and processed these samples with our GEANT-based detector simulation. As we 
analyze each reconstructed event, we reweight the event to correspond to a particular cal- 
culation for the form factors involved in the decay. This procedure allowed us to sample a 
variety of form factor calculations. For each mode, we determine the efficiency matrix for 
reconstructed versus true q 2 . Given our resolution and binning, the matrix is essentially 
diagonal, as Table II shows for the %lv form-factor calculation of Ball and Zwicky (Ball'01) 
[27]. 

For these results, we have examined the following form factors for the signal modes 
and cross-feed rates. For ittv. Ball and Zwicky (light-cone sum rules) [27], ISGW II (a 
nonrelativistic quark model) [31], and the skewed parton distributions (SPD) of Feldmann 
and Kroll [39]. Other LQCD and LCSR calculations are also considered in extracting \V u b\. 
For piv: Ball and Braun (light-cone sum rules - Ball'98) [20], ISGW II, Melikhov and Stech 
(a relativistic quark model - Melikhov'00) [38], and UKQCD (a LQCD calculation - UKQCD 
98) [8]. For r/iu, we have only considered the ISGW II form factor. The above choices for 
ixiv and piv bracket the extremes in the variation of the shape of dT / dq 2 and hence provide 
a conservative estimate of the theoretical uncertainty on the branching fractions. In general, 
the theory references provide minimal guidance on the theoretical uncertainty in the form- 
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TABLE II: The efficiency matrix in percent describing the probability that an event from a given 
generated q 2 interval reconstructs in a given q 2 interval for B° — ► tt~1 + i/ events that pass all 
cuts and reconstruct within the "signal region" of AE versus M m i u . The efficiencies are based on 
Ball'01. 



true q 2 


reconstructed q 2 


(GeV 2 ) 


- 8 8 - 16 > 16 


0-8 


2.5 0.07 0.001 


8-16 


0.07 4.6 0.06 


> 16 


0.000 0.15 4.4 



factor shapes, and the variation among the chosen calculations appears larger than the 
variation expected within a given calculation. For nominal yields and figures, we use Ball'01 
for the 7r modes and Ball'98 for the vector modes. 

We fit all the signal modes simultaneously. The parameters for the three tt~£ + u q 2 
intervals, the three p~i + v q 2 intervals, and the total r\iv branching fraction floated as free 
parameters in the fit, for a total of 7 signal parameters. The isospin and quark symmetry 
relations T(B° -> tj-^v) = 2Y(B + -> and T(B° -> p~l+v) = 2Y(B + -> p°£ + u) = 

2Y(B + — > luI + v) constrain the rates for B + relative to B°, and are assumed to hold for 
each q 2 region. We combined the three ujtv rate predictions that result from the quark 
symmetry assumption and the three piv rates to obtain the fit prediction for the total 
observed reconstructed ujIv yield. As mentioned above, only this integrated yield for uolv 
contributes to the likelihood. The two rj submodes are tied to the total r]£ + v branching 
fraction by the measured rj branching fractions and the submode reconstruction efficiencies. 
To implement the isospin constraints, we assume equal charged and neutral B production, 
/ + _ = /oo, and input a lifetime ratio of 1.083±0.017 [58]. For self-consistency, the cross-feed 
rates are constrained to the observed yields. 

The b — > c normalization in the fit varies independently for each mode, and within each 
mode for AQ = and \AQ\ = 1. The normalizations obtained are generally within 10% of 
those derived from luminosity and cross sections. The nominal fit therefore has an additional 
11 free parameters for these normalizations. 

We float the overall normalization of the generic B — > X u lv feed-down sample, deter- 
mining it from the fit. To help in determining that normalization, we take advantage of 
CLEO's recent measurement [57] of the branching fraction for b — > uiv decays with lep- 
tons in the 2.2 - 2.6 GeV/c momentum range: B(B -> XJu,2.2 < P e < 2.6 GeV/c) = 
(2.30 ± 0.38) x 10~ 4 (the "end-point branching fraction"). We constrained the B — > X u lv 
feed-down normalization by adding a x 2 term to the log likelihood of the fit: 

-21n£ - -21n£+ ^ ~ (5) 

where £> em is the measured end-point branching fraction, a cm is the total experimental uncer- 
tainty on that measurement and B cp is the branching fraction implied by the fit parameters. 
The fit prediction in each iteration is given by 

N q2 (m) 

B cp = B u n v f u n v + ^2 B mt if m ,i, (6) 

m i=l 
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where m G (7r + , 7r°, p + , p°, u, rj), B m ^ is the branching fraction for the decay mode m and 
q 2 interval % in that iteration, f m ^ is the fraction of charged leptons for that mode and q 2 
interval that are predicted by the form-factor calculation to lie in the end-point region, E u t v 
is the branching fraction for the B — > X u tv feed down background in that iteration, and 
fu£ u is the fraction of charged leptons in the end-point momentum range obtained from our 
model. 

The systematic error evaluation for the B — > X u tv feed down, and checks using alterna- 
tive procedures, are described below. The normalization is floated independently for each 
systematic variation of the various Monte Carlo, continuum, or fake samples described below 
so that the effect on the background normalization of mismodeling within the simulation is 
properly assessed. 

In summary, we have nineteen free parameters in the fit: the seven signal rates, the 
eleven generic b c background normalizations, and the one generic B — > X u tv feed- 
down background normalization. The continuum background and fake-lepton background 
samples are absolutely normalized and their rates do not float in the fit. In fits discussed 
below for which we include the \AQ\ = 1 information in the vector- meson modes, there 
are an additional 3 b — > c background normalization parameters, for a total of 22 free 
parameters. 

C. Checks and results 

We have examined the reliability of our fitting procedure via a bootstrap technique. We 
created 100 mock data samples by randomly choosing a subset of events from each of our 
Monte Carlo samples. From fits to these samples we found that our procedure reproduces 
the branching fractions without bias, and that the scatter of central values agrees with 
the uncertainties reported by the fit to better than 15%. These studies were done with 
the \AQ\ = 1 data included in the vector modes as well as in the pseudoscalar modes. 
The distribution of likelihoods that we obtained is shown in Figure 2. At the time of the 
study, we the \AQ\ = 1 data included in the vector modes. For comparison, the likelihood 
obtained from a comparable fit to the data is also shown. As discussed above, this fit has 
406 — 22 = 384 degrees of freedom. The result from the fit to data is reasonable. 

For the actual nominal fit to the data (no \AQ\ = 1 data in the vector modes), we 
obtained a value — 2 InL = 240.3 for 259 — 19 degrees of freedom. Most bins in the data fit 
have sizable statistics, so interpretation of — 21nL as a x 2 is reasonable. The probability of 
X 2 for the fit to the data is 0.48. 

In Figures 3 through 6 we show the M m £ U (AE) distributions in the AE (M m& ,) signal 
band for the individual q 2 regions examined for ixtv and for plv. For p£u, we show both 
the distributions with the nominal 1.5 GeV/c minimum lepton momentum requirement and 
with the more restrictive 2.0 GeV/c requirement of the original CLEO analysis. The fits 
describe the data in these regions well. The distributions summed over q 2 for the 7r and p 
modes and for ujIv and r\iv are shown in Fig. 7. The loiv mode remains consistent both with 
the level expected given the ptv rate and with pure background. Unless otherwise specified, 
the normalizations in all figures derive from the fit with the requirement pe > 1.5 GeV/c in 
the vector modes. 

The lepton momentum spectra and cos9 W £ distributions in the (M m& ,, AE) signal bin 
are shown in Figures 9 and 10. This information is not used in the fit, but shows good 
agreement with the signals preferred in the fit. The nn mass distribution for the combined 
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300 340 380 420 460 

-2 ln^ 

FIG. 2: Distribution of — 21nL from the bootstrap procedure described in the text. The arrow 
indicates the value obtained from the corresponding fit to the data. 

plv modes is shown in Fig. 8. 

The branching fractions from the nominal fit are summarized in Table III. The results are 
remarkably stable as the lepton momentum requirement in the vector modes is varied. The 
greatest variation is observed in the lowest q 2 interval in the piv modes, which we expected 
because of the larger role that interference between the form factors plays in that region. 

Use of a x 2 ~based fitting procedure produced similar results, though we saw clearly that 
low statistics bins had an undue influence on the results of that fitter. Such sensitivity was 
eliminated with the log likelihood minimization. 

The increase in —2 ln£ from best fit to B(B + — > r}£ + is) = is 10.4, corresponding roughly 
to a 3.2a statistical significance. 

V. EXPERIMENTAL SYSTEMATICS 

Table IV summarizes the contributions to the systematic errors for the nominal analysis. 
The dominant contribution is from uncertainties in 'V simulation," which includes inaccu- 
racies in detector simulation and uncertainty in the decay model of the nonsignal B. The 
breakdown of 'V simulation" into its component parts is given in Table V (and with lepton 
momentum cuts for vector modes of 1.75 GeV/c and 2.0 GeV/c, in Tables VI and VII, 
respectively). 
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FIG. 3: M mlv (left) and AE (right) in the AE and M m£u "signal" band requiring AQ = for the 
combined ir^jir modes. The points are the on-resonance data. The histogram components, from 
bottom to top, are b —* c (fine 45° hatch), continuum (grey or green cross hatch), fake leptons 
(cyan or dark grey), feed down from other B — » X u £u modes (yellow or light grey), cross feed 
from the vector and ij modes into the reconstructed modes (red or black fine 135° hatch), cross 
feed among the ir modes (coarse 135° hatch), and signal (open). The normalizations are from the 
nominal fit. 
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FIG. 4: Mmiu (left) and AE (right) in the AE and M mt , v "signal" band requiring |AQ| = 1 for 
the combined tt^ 1 ,^ modes. The points are the on-resonance data. See Fig. 3 for component and 
normalization descriptions. 
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FIG. 5: M mtv (left) and AE (right) in the AE and M m g v "signal" band requiring AQ = for 
the combined p^,p° modes with the requirement pe > 1.5 GeV/c in the vector modes. The points 
are the on-resonance data. The hatching and normalization are as in Fig. 3 except that the red 
or black fine 135° hatch cross feed component is from tt and rj modes into the p modes, and the 
coarse hatch cross feed component is from among the vector modes. 
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FIG. 6: M mlv (left) and AE (right) in the AE and M miu "signal" band requiring AQ = for the 
combined p^,p® modes with the requirement pi > 2.0 GeV/c in the vector modes. The points are 
the on-resonance data. The hatching is as in Fig. 5. The normalizations come from the fit with 
the corresponding lepton momentum requirement. 
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FIG. 7: M m £ u (left), AE (right) in the AE and M m i v signal bands for AQ = and summed over 
the entire q 2 range for the combined tt modes (top), p modes (row 2), uj (row 3), and r] (bottom) 
modes. See Figs. 3 and 5 for component and normalization descriptions. For rj there is only a 
single cross-feed component from the non-77 modes (red or black fine 135° hatch. 
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FIG. 8: Reconstructed mass distributions for p — > itit in the (M m £ u ,AE) signal bin for the two 
analyses with pi > 1.5 GeV/c in the vector modes (left) and with p£ > 2.0 GeV/c (right). See 
Fig. 5 for component and normalization descriptions. 
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FIG. 9: Charged lepton momentum spectrum (left) and cos Qy/i (right) for the combined iriu modes 
in the three q 2 intervals. See Fig. 3 for component and normalization descriptions. 
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FIG. 10: Charged lepton momentum spectrum (left) and cos 9wi distributions for the combined 
ptv modes in the three q 2 intervals. The cos Qy/l distribution is shown for both the analysis with 
Pi > 1.5 GeV/c in the vector modes (center) and for p? > 2.0 GeV/c (right). See Fig. 5 for 
component and normalization descriptions. 



We investigated the systematic uncertainties in V simulation" by modifying, for each sys- 
tematic contribution under consideration, the reconstruction output of all of the Monte Carlo 
samples used in the fit. Using independent studies by CLEO for this and other analyses, 
our modifications reflected the uncertainties in charged-particle-finding and photon-finding 
efficiencies, simulation of false charged particles and photons, charged particle momentum 
resolution, photon energy resolution, hadronic shower simulation, and charged particle iden- 
tification. In addition, we reweighted the Monte Carlo samples to account for the uncertain- 
ties in the rate and spectrum for K\ production in B decay and in the process b — > c — > slu, 
both of which affect the background rate into the signal region. The full MC samples were 
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TABLE III: Summary of branching fractions from the nominal fit using the Ball'01 and Ball'98 
form factors for the ir and p modes, respectively. The first uncertainties are statistical and the 
second systematic (see Section V). The results for the fits with more restrictive lepton momentum 
requirements in the vector modes are also shown. The q 2 intervals are specified in GeV 2 . 

Mode B n 2 

interval analysis requirement (vector modes) 

xlO 4 pi>1.5GeV/c pi>1.75GeV/c p e > 2.0 GeV/c 

B° n-£ + v iBtotai 1-33 ± 0.18 ± 0.11 1.31 ± 0.18 ± 0.11 1.32 ± 0.18 ± 0.12 

B <8 0.43 ± 0.11 ± 0.05 0.43 ± 0.11 ± 0.05 0.42 ± 0.11 ± 0.05 

B 8 - W 0.65 ± 0.11 ± 0.07 0.65 ± 0.11 ± 0.07 0.66 ± 0.11 ± 0.07 

B> 16 0.25 ± 0.09 ± 0.04 0.24 ± 0.09 ± 0.04 0.24 ± 0.09 ± 0.05 

B° -+ p-l+v B totSkl 2.17 ± 0.34 ±£54 2.34 ± 0.34 ±£51 2.29 ± 0.35 ±£49 

B <8 0.43 ± 0.20 0.50 ± 0.20 t° Q f 2 0.62 ± 0.22 ±g;|§ 

B 8 - W 1.24 ± 0.26 to 2 3 7 3 1.32 ± 0.26 1.11 ± 0.25 ±g;|| 

g>i 6 0.50 ± 0.10 ±g;gf 0.52 ± 0.10 ±g;gg 0.56 ± 0.10 Ijjjg 

B+ -> r]£+u B totSL i 0.84 ± 0.31 ± 0.16 0.84 ± 0.31 ± 0.16 0.83 ± 0.31 ± 0.15 



TABLE IV: Contributions to the systematic error (%) in each total and partial branching fraction 
(B). Simulation of the detector and the second B contribute to v simulation. 

■ntv p{u)lv 
q 2 interval (GeV 2 ) q 2 interval (GeV 2 ) 



Systematic 


#total 


< 8 


8-16 


> 16 


#total 


< 8 


8-16 


> 16 


7] 


v simulation 


6.8 


10.5 


9.2 


17.2 


18.7 


41.7 


19.4 


13.5 


17.3 


B -► D/D*/D**/D NR X£v 


1.7 


2.5 


1.9 


3.2 


2.0 


21.4 


4.7 


4.2 


5.5 


B — > X u tv feed down 


0.5 


3.0 


1.8 


1.9 


8.3 


23.8 


6.1 


5.6 


1.6 


Continuum smoothing 


1.0 


2.0 


0.2 


2.0 


3.0 


10.0 


1.0 


2.0 


2.0 


Fakes 


3.0 


3.0 


3.0 


3.0 


3.0 


3.0 


3.0 


3.0 


3.0 


Lepton ID 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


/+- //oo 


2.4 


2.6 


2.3 


2.2 


0.0 


2.5 


1.0 


0.1 


4.1 


T B+/ T B" 


0.2 


0.1 


0.3 


0.5 


2.1 


4.2 


1.4 


2.1 


1.4 


Isospin 


0.0 


0.0 


0.0 


0.2 


2.4 


1.9 


2.7 


2.3 


0.1 


Luminosity 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


Upper 


8.6 


12.4 


10.7 


18.3 


21.4 


53.9 


21.5 


16.2 


19.3 


Non Resonant 










-13 


-9 


-15 


-14 




Lower 


8.6 


12.4 


10.7 


18.3 


25.1 


54.7 


26.2 


21.4 


19.3 



re-analyzed for each variation to allow for leakage of events across the selection boundaries. 
The variations are described in more detail in Appendix A. 

For many of the variations in the simulation, we expect a cancellation between the change 
in the signal yield and the change in the efficiency. (Note that we are not changing the 
analysis - the data yields remain unchanged.) The cancellation arises as follows. If we 
degrade the reconstructed neutrino, the efficiency for signal is reduced, but background 
tends to smear more readily into the signal region. Hence the signal yield also tends to be 
reduced, offsetting the change in efficiency. Because of the expected imperfections in our 
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simulation, we do not expect the observed cancellation to be perfectly reliable. For each 
variation, we therefore assign an additional uncertainty in the branching fraction so that the 
total fractional uncertainty estimate is 

o" = o"br © ^ min(a yield , a eS ). (7) 

In this expression, <tbr is the percentage change in the branching fraction from the fit, Oy ic i d 
is the percentage change in the "signal bin" yield, and <T e ff is the percentage change in the 
"signal bin" efficiency For complete cancellation (cr yic i d = a eS ; ct B r = 0), the additional 
term amounts to the addition in quadrature of one third of the change observed in the yield 
and in the efficiency When no cancellation is expected, the additional term is zero. The 
values for <r y i C id and <7 e fr are estimated by examining the changes in the "signal bin." 

Note that because of correlations between the three q 2 intervals in a given mode, the 
sum of the modes tends to be less sensitive to the systematic variations than the individual 
intervals themselves. 

Consider now the items in Table IV other than 'V simulation." We reweight the Monte 
Carlo to allow variation in the relative rates for D£u, D*£u, and (Dnirjiv, both for res- 
onant Dnn and nonresonant Dnn. We vary the rates by ±8%, ±6%, ±30%, and ±30%, 
respectively. Note that if we completely eliminate any one of these charmed modes except 
D*lv, the total branching fractions for n and p remain stable within 4% of themselves, which 
demonstrates that we are quite insensitive to the details of the poorly measured nonresonant 
and resonant (Dnn) modes. Zeroing D*iv completely causes changes of only 15%, further 
demonstrating our insensitivity to the detailed modeling of the b — > civ process. 

For the B — > X u tv background, we evaluate two contributions to the systematic uncer- 
tainty. First, we vary the nonperturbative parameters of the inclusive spectrum used to drive 
the X u iv simulation within the uncertainties obtained from the B — > X s 7 analysis that were 
used in the recent end-point measurement [56, 57]. That analysis provides an error ellipse for 
the HQET parameters Ai versus A, and we choose the points on that ellipse that make the 
maximal change. The second contribution regards uncertainty in the hadronization of the 
final state light quarks. We change from our model that marries the ISGW II exclusive and 
OPE inclusive calculations (see previous section) to a purely "nonresonant" hadronization 
procedure (similar to that of JETSET [60]). The hadronization is nonresonant in the sense 
that single hadron final states (e.g., a\tv) are not produced. Resonances can appear in the 
multihadron final state (e.g., pnlv). To avoid overlap of the nonresonant sample with the 
signal modes, we eliminate B — > X u lv events with a low mass nn final state. The uncer- 
tainties presented correspond to a minimum M nn of 1 GeV. Variation of that threshold over 
the 0.9 - 1.1 GeV range results in similar systematic estimates. As a crosscheck, we have 
also used the strictly resonant description of ISGW II, which yields results consistent with 
our uncertainty estimates. 

We have used different normalization schemes for the B — > X u iv background to check 
the sensitivity of the results under the normalization procedure. If we drop the end-point 
branching fraction constraint but still allow the normalization to float, we see only minor 
shifts in the results and the end-point branching fraction predicted by the fit is within one 
standard deviation of the measured value. We have also used an iterative procedure, where 
we fix the B — > X u lv normalization in the fit, but update that normalization until the 
fit's predicted end-point branching fraction converges to the central value (and then to ±1 
standard deviation) of the CLEO measurement. This procedure also gave consistent results. 
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TABLE V: Percentage change in results for a fit with a modified simulation relative to a fit to 
the nominal MC simulation for each of the variations contributing to the simulation systematic 
uncertainty. The vector modes were analyzed with the requirement pe > 1.5 GeV/c. The last row 
shows the quadrature sum of the changes. 



variation 












p't+v 




T]lv 


total 


9 or 

q l < 8 8 


: < <r < 16 


q 2 > 16 total 


q z < 8 


8 < q 2 < 16 q 2 > 16 total 


7 eff. 


2.6 


7.0 


2.7 


9.1 


11.1 


11.9 


11.1 


10.6 


5.7 


7 resol. 


4.1 


2.9 


5.4 


2.3 


2.9 


3.7 


2.3 


4.2 


9.6 


Kl shower 


1.3 


1.0 


1.4 


1.4 


6.0 


8.4 


7.2 


1.6 


2.7 


particle ID 


1.9 


2.5 


3.0 


6.3 


8.2 


27.5 


6.9 


1.1 


0.2 


split-off rejection 


1.5 


2.9 


3.0 


5.0 


1.2 


9.4 


1.8 


2.5 


5.5 


track eff. 


3.7 


4.5 


4.2 


2.6 


8.6 


13.3 


9.5 


3.4 


9.5 


track resol. 


1.0 


1.8 


2.4 


11.2 


6.2 


12.7 


6.0 


2.7 


0.9 


split-off sim. 


0.4 


1.4 


0.5 


2.3 


1.0 


10.4 


1.0 


4.7 


6.0 


K B production 


0.2 


0.1 


0.2 


0.4 


0.1 


0.8 


0.1 


0.3 


0.1 


v production 


0.5 


3.5 


2.2 


2.0 


0.6 


15.1 


4.1 


0.9 


2.9 


Total 


6.8 


10.4 


9.2 


17.2 


18.7 


41.7 


19.4 


13.5 


17.3 



As Table IV shows, uncertainty in the B — > X u £u feed down contributes little to the 
systematic error on ixtv and r)lv. For the ptv rate, however, the contribution is substantial. 

Our nominal fit assumed equal production of charged and neutral B mesons: /+_//oo = 1- 
We varied this fraction over the one standard deviation range indicated by the recent CLEO 
result / + _//oo = 1-04 ± 0.08 [61]. The relationship enters both in the fit to implement the 
isospin constraint and in the branching fraction calculation to calculate the number of B° 
mesons. We used the measured ratio of B meson lifetimes, t b +/t b o = 1.083 ± 0.017, which 
we varied by one standard deviation to assess the associated uncertainty. The ratio comes 
into the normalization of the neutral modes versus the charged modes. We have also varied 
the isospin assumption. In the nominal fit we used a ratio of 2. For the systematic estimate 
we lowered the p + : p° ratio down to 1.43, as suggested by Diaz-Cruz [62]. The deviation 

arises from p° uo mixing coupled with the large p° width. Because of the small 77 and 

uj widths, we expect negligible deviation from the the ideal factor of two for the other two 
ratios used. 

The uncertainties related to lepton identification are estimated by varying the measured 
hadronic fake rates within their uncertainties and by applying the uncertainty in the mea- 
surement of the average lepton identification efficiency. Lepton-fake uncertainties are mea- 
sured in the data as a function of momentum using cleanly tagged hadronic samples, includ- 
ing K s -> 7i + 7i- and D*± -> tt+D , D° -> K±tx+ . 

Finally, we assessed our smoothing technique for the continuum data sample. Recall that 
we use the off-resonance data distrubution with relaxed continuum-suppression combined 
with the expected shape change over the fitted AE and M m g v region that is induced by the 
relaxation. We biased the Monte Carlo prediction for the shape change by the statistical 
uncertainty in the parameterization for each of the fitted (AE, M mlv ) distributions. The 
uncertainties come from fluctuating all distributions coherently to induce the maximum 
change. 

In the p modes, there is an additional uncertainty from the unknown contribution of non- 
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TABLE VI: Percentage change in results for a fit with a modified simulation relative to a fit to 
the nominal MC simulation for each of the variations contributing to the simulation systematic 
uncertainty. The vector modes were analyzed with the requirement pi > 1.75 GeV/c. The last row 
shows the quadrature sum of the changes. 



variation 










p-i+u 




T]lv 


total 


9 r\ r- 

q l < 8 8 


■ < q l < 16 


q 2 > 16 total 


9 no 

q l < 8 8 


< q 2 < 16 q 2 > 16 total 


7 eff. 


2.6 


6.8 


2.8 


9.3 


9.7 


8.9 


10.3 


9.0 


5.9 


7 resol. 


4.0 


2.7 


5.4 


2.4 


3.2 


4.6 


2.7 


4.1 


9.7 


Kl shower 


1.4 


1.0 


1.3 


1.7 


4.6 


4.8 


6.1 


0.5 


2.6 


particle ID 


1.8 


2.7 


3.0 


6.4 


7.8 


24.2 


6.9 


1.0 


0.0 


split-off rejection 


1.5 


2.5 


3.1 


4.7 


0.5 


1.7 


0.9 


2.4 


5.0 


track eff. 


3.7 


4.3 


4.2 


2.6 


8.4 


11.9 


9.7 


3.4 


9.7 


track resol. 


1.0 


1.8 


2.6 


11.4 


4.6 


8.1 


4.9 


1.8 


0.8 


split-off sim. 


0.4 


1.5 


0.5 


2.4 


1.1 


1.5 


0.3 


5.3 


5.3 


Kl production 


0.2 


0.1 


0.1 


0.4 


0.1 


0.7 


0.1 


0.3 


0.0 


v production 


0.5 


3.5 


2.3 


2.2 


0.8 


13.3 


3.1 


0.6 


2.7 


Total 


6.7 


10.2 


9.3 


17.4 


16.7 


33.1 


18.0 


12.2 


17.0 



TABLE VII: Percentage change in results for a fit with a modified simulation relative to a fit to 
the nominal MC simulation for each of the variations contributing to the simulation systematic 
uncertainty. The vector modes were analyzed with the requirement pe > 2.0 GeV/c. The last row 
shows the quadrature sum of the changes. 



variation 










p-l + v 




r\lv 


total 


q 2 < 8 8 


. < q 2 < 16 


q 2 > 16 total 


q 2 < 8 8 


< q 2 < 16 


q 2 > 16 total 


7 eff. 


2.6 


6.8 


2.7 


8.8 


12.3 


12.3 


14.6 


8.3 


5.9 


7 resol. 


4.2 


2.7 


5.4 


3.9 


1.3 


1.1 


2.3 


1.0 


9.3 


Kl shower 


1.3 


0.9 


1.7 


1.7 


2.4 


1.8 


3.2 


1.2 


2.6 


particle ID 


1.9 


2.7 


3.1 


6.4 


7.0 


15.6 


8.1 


1.1 


0.4 


split-off rejection 


1.7 


2.7 


3.0 


5.9 


1.8 


1.5 


2.9 


0.7 


5.7 


track eff. 


3.9 


4.3 


4.5 


2.4 


4.1 


4.7 


6.5 


1.8 


9.2 


track resol. 


1.0 


1.5 


3.0 


11.8 


3.6 


8.2 


2.4 


2.7 


1.0 


split-off sim. 


0.4 


1.6 


0.5 


3.1 


1.9 


6.6 


2.8 


3.0 


5.2 


Kl production 


0.2 


0.2 


0.1 


0.4 


0.2 


0.5 


0.2 


0.4 


0.1 


v production 


0.6 


3.5 


2.4 


2.6 


0.7 


6.3 


1.5 


0.6 


2.1 


Total 


7.0 


10.2 


9.7 


18.2 


15.7 


23.9 


18.9 


9.7 


16.7 



resonant Txniv decays. While little is known about these decays, we can provide a framework 
for limiting those contributions through the study of reconstructed 7i°7i°£ ± i' decays and the 
consideration of Bose symmetry, isospin, and angular momentum. The B — > X u tv decay 
results, before hadronization, in two final-state light quarks. These can have either isospin 
I — or I — 1. Because final-state interactions preserve isospin, a final nn state is also 
restricted to / = or I — 1. From Bose symmetry considerations, the Tin state must have 
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TABLE VIII: Comparisons of the p l + v branching fractions when the tv mode and component 
are added. The parameter a that normalizes the ifiifitv component is described in the text. The 
percentage change relative to the standard fits in Table III are indicated in parentheses below the 
branching fractions. 



analysis 


a 


B(B ->■ ptv) B q2<8 GcV 2 B 8<q2<w GcV 2 B q2>16 GeV 2 


x 2 /dof 






(KT 4 ) (10~ 4 ) (10- 4 ) 





p/>1.5GeV/c 0.25 ±0.21 1.88 ± 0.35 0.393 ± 0.209 1.055 ±0.261 0.428 ± 0.099 273.7 / (280-21) 

(-13%) (-9%) (-15%) (-14%) 

p e > 1.75 GeV/c 0.22 ± 0.18 2.06 ± 0.35 0.455 ± 0.216 1.153 ±0.259 0.455 ± 0.099 271.6 / (280-21) 

(-12%) (-8%) (-9%) (-11%) 

Pl> 2.0 GeV/c 0.18 ±0.13 2.17 ±0.36 0.669 ± 0.249 1.009 ±0.241 0.496 ± 0.097 281.1 / (280-21) 

(-5%) (8%) (-9%) (-11%) 



angular momentum L even for 1 = and L odd for 1=1. Isospin considerations then imply 

I = 1, L odd t^tt : 7T+7T- : 7r°7r° = 2:1:0 
I = 0, L even 7r ± 7r° : u+ij- : 7r°7r° = 0:2:1. 

Assuming that the L = 3,5, ... configurations are suppressed relative to the L — 1 configu- 
ration, we can use e + e" scattering data and r decay data to conclude that the I — 1, L odd 
component is completely dominated by the p. A significant nonresonant contribution would 
therefore come via the I = 0, L even channel. With the 1 = rate parameterized by a, we 
expect partial widths in the ratios 

t^tt : 7r + 7T~ : ttV = 2 : 1 + 2a : a. 

To estimate the systematic due to an unknown nonresonant ti-rIv contribution, we look for 
a component, after event selection, that could mimic a piv. To constrain such a contribution, 
we add the mode 7r 7r £z/ to the fit. Procedurally, we generate ir°ir £v using the p lineshape 
and the ptv form factors. We then perform fits with the usual isospin constraint on the 
partial widths (p ± : p° — 2 : 1) replaced with the hit ratios given above. While the most 
relevant fit for the extraction of a systematic uncertainty number has the parameter a 
floating, we also fix a = to test the fit quality under the assumption that observed 7r°7r°£z/ 
yields are consistent with cross feed from other modes and the other standard backgrounds. 

In the fits, the 7r°7r°fz/ mode is treated like the u mode. Only the sum of the three q 2 
intervals contributes to the likelihood, but the signal Monte Carlo is scaled in each q 2 interval 
separately to maintain the above 7T7T ratios from one interval to the next. Figure 11 shows the 
projection onto the m n o n o distribution for fits with and without a tt 71 £u signal component. 
Note that the fit included data only from the three bins in the range 0.485 < m^o^o < 1.055 
GeV. The fit quality is excellent when the reconstructed 7r°7r°£z/ mode is included but the 
7r°7r°£z/ signal is forced to zero. Table VIII summarizes the observed changes in the p~ 
branching fraction when we float the ir ir lis signal component. The resulting 7r°7r°£z/ yield 
is consistent with zero. The shifts in the various piv branching fractions are larger effects 
than the increase in their errors due to correlations with the 7t 7t £za We thus take the shifts 
as the estimate of the uncertainty. The pseudoscalar modes shift negligibly. 

In addition to the variations above, we have performed numerous systematic checks, 
including variation of the selection criteria and investigation of electron and muon samples 
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0.5 1.0 1.5 0.5 1.0 1.5 

M^O (GeV) 

FIG. 11: The 7r°7r° mass distribution from the reconstructed ^^i^u "signal bin" from the nominal 
fit (left) and from the fit including a ^^i^u signal component (right) as described in the text. 
The points are the on-resonance data. The histogram components, from bottom to top, are b — > c 
(fine 45° hatch), continuum (grey or green cross hatch), fake leptons (cyan or dark grey), feed 
down from other B — > X u tv modes (yellow or light grey), cross feed from the signal modes into 
the reconstructed modes (red or black fine 135° hatch), and signal (open). The normalizations are 
from the corresponding fits. 



separately. We have also investigated tighter momentum requirements in the pseudoscalar 
modes. The observed variations were in general consistent within the uncertainties resulting 
from the statistical changes. 

VI. DEPENDENCE OF BRANCHING FRACTIONS ON FORM FACTORS 

In the original measurement of the exclusive charmless branching fractions [3] , there were 
two roughly comparable contributions to the branching fraction errors from the form-factor 
uncertainties. The first contribution resulted because the efficiency varied as a function of q 2 
(inescapable with a lepton momentum cut), and the data were lumped into a single q 2 bin. 
Because we now extract the rates independently in three separate q 2 ranges, this analysis 
should see a significant reduction in this effect. The second contribution resulted because 
there was significant q 2 dependence to the cross-feed rates between the pseudoscalar and 
the vector modes. Again, since we extract the rates independently as a function of q 2 , this 
dependence should be reduced. 

We have estimated the model dependence based on changes of the branching fractions 
under variation of the form-factor calculation. The previous analysis [3] found that the error 
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TABLE IX: Branching fractions B(B° — ► tt~£ + v) obtained under variation of the tt and pjiolv 
form-factor models. Shown are the results for the total branching fraction, the partial branching 
fraction in each q 2 bin, and the — 21nL for the fit. Branching fraction uncertainties are statistical 
only. The estimated model dependence is indicated after each set of variations. All branching 
fractions are in units of 10~ 4 . The tt model variations are all presented for the analysis with the 
pe > 1.5 GeV/c requirement on the vector modes. 







I 2 


interval (GeV' 


2 ) 




tt Model p Model 


^total 


B<8 


^8-16 


£>16 


—2 In I 


Ball'01 Ball'98 


1.327 ± 0.177 0.431 ± 0.106 0.651 ± 0.105 0.245 ± 0.094 


240.3 


ISGW2 Ball'98 


1.327 ± 0.176 0.431 ± 0.107 0.660 ± 0.106 


0.236 ± 0.09 


240.7 


SPD Ball'98 


1.315 ± 0.173 0.436 ± 0.106 0.650 ± 0.105 0.229 ± 0.088 


239.8 


1.7 x RMS,,- ff 


0.01 


0.004 


0.01 


0.01 




Ball'01 Ball'98 


1.33 ±0.18 


0.43 ±0.11 


0.65 ±0.11 


0.25 ± 0.09 


240.3 


Ball'01 ISGW2 


1.41 ±0.18 


0.45 ±0.11 


0.69 ±0.10 


0.27 ±0.09 


239.4 


Ball'01 Melikhov'00 


1.30 ±0.18 


0.43 ±0.11 


0.65 ±0.11 


0.22 ±0.09 


240.2 


Ball'01 UKQCD'98 


1.36 ±0.18 


0.44 ±0.11 


0.66 ±0.11 


0.26 ±0.09 


239.3 


1.7 x RMSpFF 


0.07 


0.01 


0.03 


0.03 





on the branching fraction obtained from comparison of models was larger than that obtained 
by variation of a particular form-factor parameterization within the published uncertainties 
(when given). Tables IX and X show the variation in B(B° — > tt~£ + i>) and B(B° — > p~£ + u), 
respectively, as the tt and vector form factors are varied. We have included in the set of 
models those which have the most extreme variations in shape of dT/dq 2 . For tt£i>, we find 
that our method results in almost no sensitivity to the form factor used for the signal mode 
efficiencies. We find a larger sensitivity to the variation of the vector mode form factors 
because of cross feed from those modes. For p£u, there is almost no sensitivity to the tt£v 
form factors, but significant sensitivity to the p£v form factors. 

To assign uncertainties, we use an empirical observation from the original analysis [3]. For 
that analysis, for any given model, we varied the internal parameters to determine an error 
on the rates extracted within that model. We then denned a range of potential branching 
fractions by taking the model with the lowest result and subtracting one standard deviation 
from the variations within that model, and taking the model with the highest result and 
adding one standard deviation. Our assigned uncertainty covered 70% of this range. (Note 
that this procedure gave us a more conservative range than taking one half the spread 
among the central value of the models.) Empirically, we found that this procedure agreed 
with taking 1.7 times the RMS spread among models for all quantities that we examined. 
For these results, we therefore apply this latter procedure. The results are also summarized 
in Tables IX and X. 

For purposes of direct comparison, had we adopted the procedure used in recent p£v 
analyses by the BABAR Collaboration [63] and by CLEO 2000 [4], we would assign (abso- 
lute) uncertainties of 0.06 x 10 -4 (rather than 0.07 x 10~ 4 ) and 0.33 x 10 -4 (rather than 
0.41 x 10~ 4 ) for the p~ £ + v form-factor dependence on the total branching fraction for tt~£ + v 
and p~£ + v, respectively. The p~£ + u number, 0.33 x 10~ 4 , is about half of the size seen in the 
recent BABAR measurement, which, like the CLEO 2000 measurement, is mainly sensitive 
to the end-point region pe > 2.3 GeV/c. 
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TABLE X: Branching fractions B(B° — > p~£ + u) obtained under variation of the tt£u and pjiolv 
form-factor models. Shown are the results for the total branching fraction, the partial branching 
fraction in each q 2 bin, and the — 21nL for the fit. Branching fraction uncertainties are statistical 
only. The estimated model dependence is indicated after each set of variations. All branching 
fractions are in units of 10~ 4 . The tt model variations are all presented for the analysis with the 
pn > 1.5 GeV/c requirement on the vector modes. For the vector mode form-factor variation, we 
present the results for all three momentum requirements. 







Q 2 


interval (GeV 2 ) 




7r Model p Model 


£>total 




£>8-16 


£>l 6 


— ZmL 


t -) h)ai i -) n?rvo 

Ball 01 Ball 98 


2.172 ±0.338 0.429 ±0.198 


1.244 ± 0.256 0.499 ± 0.097 


240.3 


1SGW2 Ball 98 


2.176 ± 0.338 0.430 ± 0.198 1.248 ± 0.256 0.499 ± 0.098 


240.7 


brU Bail 98 


2.169 ± 0.338 0.420 ± 0.198 1.250 ± 0.256 0.499 ± 0.098 


239.8 




0.01 


0.01 


0.004 


0.004 




pe > i.o tote v /c 












Ball 01 Bali 98 


2.17 ±0.34 


0.43 ± 0.20 


1.24 ±0.26 


0.50 ±0.10 


240.3 


Ball 01 1SGW2 


1.91 ±0.28 


0.30 ±0.13 


1.14 ±0.23 


0.47 ±0.10 


239.4 


Ball 01 Mehkhov 00 


2.56 ±0.37 


0.33 ±0.15 


1.49 ±0.31 


0.75 ±0.14 


240.2 


Bail 01 UKyCU 98 


2.08 ±0.32 


0.39 ±0.17 


1.21 ±0.25 


0.49 ±0.10 


239.3 


1.7 x RMSpFF 


0.41 


0.09 


0.22 


0.19 




p e > 1.75 GeV/c 












Ball'01 Ball'98 


2.34 ±0.34 


0.50 ±0.20 


1.32 ±0.26 


0.52 ±0.10 


241.6 


Ball'01 ISGW2 


2.03 ±0.28 


0.34 ±0.13 


1.20 ±0.23 


0.49 ±0.10 


240.3 


Ball'01 Melikhov'00 


2.74 ±0.37 


0.38 ±0.16 


1.58 ±0.31 


0.78 ±0.14 


241.4 


Ball'01 UKQCD'98 


2.23 ±0.32 


0.45 ±0.18 


1.28 ±0.25 


0.51 ±0.10 


240.4 


1.7 x RMSpFF 


0.44 


0.11 


0.24 


0.20 




Pe > 2.0 GeV/c 












Ball'01 Ball'98 


2.29 ±0.35 


0.62 ±0.22 


1.11 ±0.25 


0.56 ±0.10 


244.2 


Ball'01 ISGW2 


1.89 ±0.27 


0.38 ±0.13 


0.98 ± 0.22 


0.54 ±0.09 


243.4 


Ball'01 Melikhov'00 


2.66 ±0.38 


0.48 ±0.17 


1.36 ±0.31 


0.83 ±0.14 


244.6 


Ball'01 UKQCD'98 


2.15 ±0.32 


0.54 ±0.19 


1.07 ±0.24 


0.55 ±0.09 


243.3 


1.7 x RMSpFF 


0.47 


0.15 


0.24 


0.21 





We stress that the form factors from any given model are not used to constrain the relative 
rates extracted in each of the three q 2 regions. Only the efficiencies within each q 2 range are 
modified. Hence the quality of the fit used to extract the rates does not discriminate among 
different form-factor descriptions. This discrimination is discussed in the following section. 

Overall, our procedure has drastically reduced the sensitivity of the niv result to both 
the 7r£z/ and the vector-mode form factors. There is essentially no dependence on the Tciu 
form factors themselves. The combined sensitivity to both the n and p form factors is about 
one third that of the previous CLEO nlv analysis. 

The plv variation remains significant, though again this analysis shows essentially no 
dependence on the niu form factor. The overall uncertainty of the form factors has reduced 
to about 80% of the original CLEO piv measurement [3] (which had a smaller form-factor 
dependence than the 2000 CLEO piv analysis [4]). As one tightens the lepton momentum 
requirement, the model dependence increases slightly over the range we have studied. As 
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expected, the lowest q 2 interval shows the greatest sensitivity (fractionally) to the variation 
in the range. For a given model, the variation of the total branching fraction as the lepton 
momentum requirement is varied is small compared to the variation among models for a 
given momentum requirement. (The RMS variation of the former is about 30% of the RMS 
variation of the latter.) We speculate that the dominant model dependence likely arises 
from our cos6y\/£ > requirement, which we applied to suppress b — > c background. Either 
finer q 2 binning or an alternate means of background suppression would provide a route for 
further reduction of the form-factor dependence. 

For the r/iis branching fraction, we find a dependence of 0.04 x 10" 4 from variation of the 
nlv form factors and 0.01 x 10~ 4 from variation of the plv form factors. The only r\lv form 
factor that we consider is ISGW II [31]. However, the f]£u analysis presented here is almost 
identical to the original i\iv analysis. We therefore take the form-factor dependence of 10% 
found in that analysis as an estimate of the uncertainty from the rj form factors. As the ptv 
form factors contributed substantially to the 10% uncertainty in the previous analysis, yet 
contribute negligibly to r]£v, the 10% should be a conservative estimate. 

The results presented here agree well with the previous CLEO measurements and the 
recent BABAR piv measurement. The results of the original CLEO measurement [3] are 
superseded by this measurement. The results of the CLEO 2000 measurement [4] are essen- 
tially statistically independent of those presented here. 

VII. EXTRACTION OF \V ub \ AND DISCRIMINATION OF MODELS 

We extract \V u b\ from the measured rates for n£v only, for ptv only, and then by using 
the combined information from those two modes. In all cases, the |V4&| extraction is based 
on the results from the the analysis requiring pe > 1.5 GeV/c in the vector modes. We use 
a B° lifetime of (1.542 ± 0.016) ps [58]. 

A. \Vub\ from B — ► -nlv 

For tt£v, we first explore fitting q 2 distributions from various form-factor predictions to the 
measured rates in the three q 2 bins. To be self-consistent, we extract \V u b\ for a particular 
form factor using the rates from the fit with that model. In practice, as we have seen, 
this makes little difference in the 7r modes in this analysis. Since each model predicts the 
total rate modulo \V u b\, \V u b\ becomes the one free parameter for the fit that normalizes the 
prediction to the observed rates. The quality of the fit measures how well the form-factor 
shape describes the data, so provides one means of discrimination among form factors. The 
results of this procedure are summarized in Table XI. For the three calculations that have 
been used for both efficiency and \V u b\ extraction, the data rates with the best fits for each 
predicted form factor are shown in Figure 12. The probability of \ 2 i n our various fits for 
the ISGW II model varies between one to three percent, indicating that this model is likely 
to be less reliable for determination of \V u b\ from niv. Note further that the spread among 
the central values from the various calculations is fairly small relative to the uncertainties 
quoted in the calculations themselves. 

Because the extracted rates in the q 2 intervals are now essentially independent of the n£v 
form factor, one can extract \V u b\ from our results for form factors not considered here. We 
provide in Appendix B a detailed methodology for doing so. 
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FIG. 12: Measured branching fractions in the restricted q 2 intervals for B° — > ir~l + v (points) and 
the best fit to the predicted dT/dq 2 (histograms) for the three models used to extract both rates 
and \V u b\- The data points have small horizontal offsets introduced for clarity. The last bin has 
been artificially truncated at 24 GeV 2 in the plot - the information out to q 2 ^^ has been included 
in the work. 



To determine the effect of the systematic uncertainties, we repeat the above fit using 
the three q 2 rates obtained from the branching ratio fit after each systematic variation. 
This procedure automatically accounts for correlations among the three intervals. We then 
increase the uncertainty for each variation by one half of the fractional error introduced by 
the second term in Equation 7. The factor of one half arises from the square root involved 
in extraction of \V u b\ from the rate. 

As we discuss below, each of the form-factor calculations used to extract \V u b\ from the full 
q 2 range has some measure of model dependence. We determine a systematic error in \V u b\ 
from the quoted theoretical uncertainty in form-factor normalizations, with the following 
procedure. For each form factor used, we recalculate \V u b\ when we increase or decrease 
the form-factor normalization by one standard deviation. Due to the poor agreement of the 
ISGW II form factor with the niv data in conjunction with the somewhat ad hoc assumptions 
about the form-factor g 2 -dependence in that mode, we drop ISGW II from consideration. 
From the others, we find the minimum value V m i n and the maximum value V mSuX . We then 
assign an asymmetric error of 70% of the deviation relative to the nominal central value - 
that is, we take 0.7(V A max — V nom ) and 0.7(Vn O m — Knin)- Because the result obtained using 
Ball'01 is close to the mean, we take that result as the nominal value. Note that when a 
symmetric theory error is quoted on the rate, we re-interpret that error as symmetric on 
the amplitude. To be precise, we map 7th ± o"th to 7th ± o"th(l ± 0"th/(47th))- This procedure 
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TABLE XI: \V u b\ extracted from fits to the rates measured in the three q 2 intervals for a variety 
of form factors for ir£v. The table indicates form-factor calculation, \V u b\ with statistical error 
only, predicted r^ h /|V^&| 2 with the estimated theoretical uncertainty, the x 2 f° r the fit, and the 
probability of x 2 given the two degrees of freedom. 



7r Model 


p Model 


\v uh \ x io 3 r^ h 


l\V ub ? (PS" 1 ) 


Fit x 2 


P(x 2 ) 


Ball'01 


Ball'98 


3.21 ±0.21 


8.4l 3 ;> 


1.0 


0.61 


KRWWY" 


Ball'98 


3.40 ± 0.23 


7.3 ±2.5 


5.3 


0.07 


ISGW2 


Ball'98 


2.90 ± 0.20 


9.6 ±4.8 


7.3 


0.03 


SPD 


Ball'98 


2.96 ±0.19 


9.6 ±2.9 


4.0 


0.14 


Ball'01 


Ball'98 


3.21 ±0.21 


8.4±- 


1.0 


0.61 


Ball'01 


ISGW2 


3.31 ±0.20 


8.4l|l 


1.2 


0.55 


Ball'01 


Melikhov'00 3.18 ±0.21 


8.4l 3 i 


0.9 


0.63 


Ball'01 


UKQCD'98 


3.24 ±0.20 


8.4l|l 


1.1 


0.59 



a Uses rates determined with the Ball'01 form factor, 
yields 

\V ub \ = (3.21 ±0.21 ±0.14 toll ±0.10) x 10" 3 , (8) 

where the errors are statistical, experimental systematic, the estimated uncertainties from 
the 7i£u form-factor shape and normalization, and the piv form-factor shape, respectively. 
The piv form-factor contribution has been estimated using the 1.7<trms prescription. 

Again for direct comparison with other experiments, taking one half, rather than 70%, 
as the scale factor for estimating the uncertainties yields \V u b\ = (3.21 ± 0.21 ± 0.14 I0I32 

± 

0.07) x 10" 3 . 

Note that the error on \V u b\ from the uncertainty in the rates under variation of form 
factors is completely dwarfed by the error arising from uncertainty in the theoretical nor- 
malization of the form factor. 

Our second, preferred, method for determining \V u b\ attempts to reduce the number 
of modeling assumptions and hence to provide a more robust uncertainty estimate. We 
therefore limit our consideration to form factors determined from LCSR and from LQCD 
calculations are QCD-based calculations. These calculations, however, are only valid over a 
restricted q 2 region. The LCSR assumptions are expected to break down for q 2 > 16 GeV 2 , 
while the current LQCD calculations are valid only for q 2 > 16 GeV 2 . Extrapolation outside 
of these ranges therefore introduces a dependence on the form used for the extrapolation. 
This introduces another uncertainty that is difficult to asses. To minimize this uncertainty, 
we extract \V u b\ from these more restricted regions. For LQCD, we determine \V u b\ from the 
measured rate and the calculated rate in the range q 2 > 16 GeV 2 . For LCSR, we determine 
\V u b\ by fitting the calculated LCSR rates to the measured rates in the two q 2 intervals below 
16 GeV 2 . The results are shown in Table XII. 

To produce a final LQCD result for the q 2 > 16 GeV 2 region, we take a statistically- 
weighted average of the different LQCD results. To the precision quoted, we obtained 
identical results if we based the statistical weights on the upper, the lower, or the average of 
the asymmetric statistical errors quoted in Table XII. We assume the systematic errors are 
completely correlated among the different calculations: if ctj is the statistical weight used in 
the average for calculation % and <jj is the fractional systematic error for that calculation, 
then the total fractional systematic error a assigned to the average is a = c^i^i- The 
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TABLE XII: Values for \V u b\ obtained using form factors (FF) from light-cone sum rules in the q 2 
interval — 16 GeV 2 (top two rows) and from LQCD for q 2 > 16 GeV 2 (bottom five rows). Only 
the statistical errors on \V u b\ are indicated. The data rates obtained using Ball'01 for ixlv and 
Ball'98 for ptv were used as the input for all values obtained. 



vr FF \V uh \ x 10 3 Tf t 


'\V ub \ 2 (P^ 1 ) 


Fit x 2 


P(x 2 ) 


Ball'01 3.20 ± 0.22 


6-9™ 


1.0 


0.32 


KRWWY 3.46 ± 0.24 


5.7 ± 1.9 


5.0 


0.025 



FNAL a 2.88 ±0.55 1.911^ ± 0.31 

JLQCD 6 3.05 ±0.58 1.71±g;|| ± 0.46 

APE C 2.97 ±0.57 1.80±g;|? ± 0.47 

UKQCD d 2.63 ± 0.50 2.3lg;^ ± 0.51 

average 6 2.88 ± 0.55 1.92±g;?| d= 0.47 



a The authors of [16] have provided the rate integrated over this range and the corresponding 
uncertainty. 

b The authors of [17] have provided the rate integrated over this range and the corresponding 
uncertainty 

c We have integrated over the restricted q 2 interval to obtain rates using the FF parameterization 
from the two APE methods, scaled the uncertainties accordingly and performed a simple average 
of the two rates. 

d We have integrated the FF parameterization over the restricted q 2 interval to obtain the central 
value and have scaled the uncertainties accordingly. 
e See text. 



theoretical errors quoted in Table XII do not include any uncertainty from the quenched 
approximation, which is estimated to be in the 10% to 20% range. We add an additional 15% 
in quadrature to the systematic uncertainty just described to obtain the average theoretical 
systematic uncertainty quoted in the table. 

From our average of the LQCD-based results, we estimate 

\V ub \ q *>i 6 Gc V 2 = (2.88 ± 0.55 ± 0.30 ± 0.18) x 10" 3 , (9) 

where the errors are statistical, experimental systematic, LQCD uncertainties, and ptv form- 
factor dependence, respectively. The LQCD uncertainties have been combined in quadrature. 

Taking the simple average of the two LCSR values and again using the 70% range to 
estimate the theoretical uncertainty, we characterize the LCSR results as 

|Kt|,»<i6 Gc V 2 = (3.33 ± 0.24 ± 0.15 t° H ± 0.06) x 10" 3 . (10) 

Using the fractional errors from the LCSR calculations alone gives similar theoretical un- 
certainties. 

We average the LQCD and LCSR results, with correlated experimental systematics taken 
into account, according to the procedure laid out in Appendix C. The LQCD value enters 
the average with a weight of a n = 0.20. As noted in the appendix, we choose the weight to 
minimize the total overall uncertainty. To be conservative, we have treated the theoretical 
uncertainties as if they were completely correlated. 

\V ub \ = (3.24 ± 0.22 ± 0.13 t° f 9 ± 0.09) x 10" 3 . (11) 

We take this as the more reliable determination of \ V u b\ from our complete data in this mode. 
The variations in |14&| and our averages are illustrated in Figure 13. 
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FIG. 13: Values for \V u b\ obtained from ir£v using the entire q 2 range for the various form- factor 
calculations (top block), using LQCD for q 2 > 16 GeV 2 (second block), using LCSR for q 2 < 16 
GeV 2 (third block), and our average of the last two (bottom block) for -kIv only and for -kIv and 
plv combined. In all cases, the top bar indicates the statistical and all the experimental systematics 
(combined in quadrature), the lower bar indicates the approximate "one standard deviation" range 
of motion due to the theoretical uncertainties. 

B. | Kb | from B — ► ptv 

We proceed with B — > piv in much the same fashion as with B — > nlv. The fits of 
the different form factors to the rates extracted from the three q 2 intervals in the data are 
illustrated in Figure 14 and are summarized in Table XIII. Because of the relatively large 
variation in the rates extracted from the data using the different form-factor calculations, 
we again perform the extraction of \V u b\ entirely within the context of a given form-factor 
calculation. In general, the theoretical predictions do not match the data as well as we saw 
for the -nlv mode. In spite of some of the poor fits, we consider all four sets of form factors 
as we estimate |T4fc| with this mode. As we expected from the branching fraction results, 
the \V u b\ extracted from the piv information does not depend on the ntv form factor used 
in the analysis. 

For an estimate of \V u i\ based on the models and fits in Table XIII, we take the Ball'98 
results as the central value. Estimating the uncertainties as described in the previous section, 
we obtain 

\V ub \ = (2.90 ± 0.21 t° H t° 7 i) x lO" 3 , (12) 

where the errors are statistical, experimental systematic, and the estimated uncertainties 
from 70% of the total spread in the results as we vary the piv form-factor calculations over 
±1 standard deviation, respectively. This estimate is similar to, though somewhat larger 
than, that obtained from the quoted Ball'98 uncertainty. 



33 



3600103-005 





1 R 
I .O 


■<3- 


1 .6 


o 




X 


1.4 


-t— ' 

o 




U 


1.2 


CM 

cr 


1 .0 






cr 


0.8 


35 






0.6 








0.4 




0.2 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



T --£- ■ 



% 

i : I i 



• ; MelikRov -] 

■ Ball '98 

▲ ISGW II 

O UKQCDW 



i 



ii 

A. 



..L.J ' 



i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i 







10 

q 2 (GeV'') 



15 



20 



FIG. 14: Measured branching fractions in the restricted q 2 intervals for B° — > p~i + v (points) and 
the best fit to the predicted dT/dq 2 (histograms) for the three models used to extract both rates 
and \V u b\. The data points have small horizontal offsets introduced for clarity. 



TABLE XIII: \V u b\ extracted from fits to the rates measured in the three q 2 intervals for a variety 
of form factors for plv. The table indicates form-factor calculation, \V U \\ with statistical error 
only, predicted r* h /|I4ft| 2 with the estimated theoretical uncertainty, the x 2 f° r the fit, and the 
probability of x 2 given the two degrees of freedom. 



7T Model 


p Model 


\Vub\ x 10 3 


C/lKfcl 2 (ps- 


!) Fit x 2 


P(x 2 ) 


Ball'01 


Ball'98 


2.90 ±0.21 


16.9 ±5.1 


7.6 


0.02 


BalPOl 


ISGW2 


2.96 ±0.21 


14.2 ±7.1 


3.3 


0.19 


BalPOl 


Melikhov'OO 2.46 ± 0.17 


26.2 ±5.2 


8.1 


0.02 


Ball'Ol 


UKQCD'9£ 


] 2.88 ±0.20 


16.511* 


5.2 


0.08 


BalPOl 


Ball'98 


2.90 ±0.21 


16.9 ±5.1 


7.6 


0.02 


ISGW2 


Ball'98 


2.90 ±0.21 


16.9 ±5.1 


7.6 


0.02 


SPD 


Ball'98 


2.90 ±0.21 


16.9 ±5.1 


7.8 


0.02 



Restricting ourselves to the theoretically more reliable use of LQCD for q 2 > 16 GeV 2 and 
LCSR for q 2 < 16 GeV 2 , we have only the two results listed in Table XIV. In addition to the 
theoretical uncertainty quoted for UKQCD'98, we add an additional 20% in quadrature as 
an estimate of the quenching uncertainty. This is larger than for the tt£u case both because 
the p is a broad resonance and because of the potential for larger biases from quenching 
given the interference between the various form factors. We also apply our reinterpretation 
of symmetrical theoretical errors on the rate as symmetric errors on the amplitude. The 
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TABLE XIV: Values for \V ub \ obtained using form factors (FF) from light-cone sum rules in the 
q 2 interval — 16 GeV 2 (first row) and from LQCD for q 2 > 16 GeV 2 (second row). Only the 
statistical errors are indicated. The data rates obtained using Ball'01 for ir£v and Ball'98 for plv 
were used as the input for all values obtained. 



p FF \V ub \ x 10 3 


C/lKfcl 2 (ps- 1 ) 


Fit x 2 


P(x 2 ) 


Ball'98 2.67 ± 0.27 


14.2 ±4.3 


4.5 


0.03 


UKQCD'98 3.34 ±0.32 


9 n+U.t)2 
^• y -0.40 







results in the two q 2 intervals are thus 

\V ub \ q 2> w GcV2 = (3.34 ± 0.32 1°; 2 6 7 ±g;tS) x 10~ 3 , (13) 

and 

\V ub \ q 2 <w GcV2 = (2.67 ± 0.27 ) x 10- 3 . (14) 

We average the LQCD and LCSR results, with correlated experimental systematics taken 
into account. We again employ the procedure described in Appendix C. The optimal 
weight for combining the two intervals, treating the systematic uncertainties as completely 
correlated, is a p = 0.5. 

\V ub \ = (3.00 ± 0.21 ±g;i§ t° f 8 ± 0.28) x lO" 3 . (15) 

The errors are statistical, experimental systematic, theoretical systematic based on the 
LQCD and LCSR uncertainties, and piv form-factor shape uncertainty. To be conserva- 
tive, we have assigned the latter error based on the variation seen in the total branching 
fraction in this mode. The contribution from the niv form-factor shape is negligible. Again, 
we take this as our preferred method of extracting \ V ub \ from our piv data. 
The |T4b| results obtained from ptv are shown in Figure 15. 

C \V u b\ from a Combination of B — > tv£v and B —>■ plv 

We have averaged the \V ub \ determinations obtained separately from the B — > niv and 
B — > plv modes. For this average, we considered only the results obtained using the LCSR 
and LQCD calculations applied to the q 2 < 16 GeV 2 and q 2 > 16 GeV 2 results, respectively. 
The averaging procedure amounts to the determination of the optimal weight f3 to be applied 
to the LCSR and LQCD average obtained from B — > ixlv relative to that obtained from 
B — > piv (see Appendix C). We held the values a n and a p , each of which determines the 
weight of the LQCD result relative to the LCSR result in the individual mode, fixed at the 
optimal values found in the preceding subsections. The weight (3 = 0.7 provided the optimal 
combination. With this weighting, we find 

|K&| = (3.17 ± 0.17 lg;l?lg;i± 0.03) x lO" 3 . (16) 

The errors are statistical, experimental systematic, theoretical systematic based on the 
LQCD and LCSR uncertainties, and piv form-factor shape uncertainty, respectively. Note 
that because of cross feed among the modes considered, the tt£v and piv modes are anticor- 
related, resulting, in particular, in the minimal dependence of the average result on the piv 
form-factor shape. 
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FIG. 15: Values for \V u b\ obtained from plv using the entire q 2 range for the various form- factor 
calculations (top block), using LQCD for q 2 > 16 GeV 2 (second block), using LCSR for q 2 < 16 
GeV 2 (third block), and our average of the last two (bottom block) for plv only and for nlv and piv 
combined. In all cases, the top bar indicates the statistical and all the experimental systematics 
(combined in quadrature), the lower bar indicated the approximate "one standard deviation" range 
of motion due to the theoretical uncertainties. 

VIII. SUMMARY 

With a sample of 9.7 x 10 6 BB pairs, we have studied B decays to niu, p£v, u£v, and r]£v, 
where £ = e or p. From the combination of a broad momentum range for the charged lepton 
momentum and independent extraction of rates in three separate q 2 intervals, we were able 
to reduce the uncertainties from modeling within the form-factor calculations. For the decay 
B° — > tv^£ + v, we have determined the branching fractions 

8(0 < q 2 < 8 GeV 2 ) = (0.43 ± 0.11 ± 0.05 ± 0.004 ± 0.01) x 10" 4 (17) 
B(8 < q 2 < 16 GeV 2 ) = (0.65 ± 0.11 ± 0.07 ± 0.01 ± 0.03) x 10~ 4 
B(q 2 > 16 GeV 2 ) = (0.25 ± 0.09 ± 0.04 ± 0.01 ± 0.03) x 10~ 4 . 

Combining these rates and taking into account correlated systematic uncertainties, we obtain 

B(B° -> ix- tv) = (1.33 ± 0.18 ± 0.11 ± 0.01 ± 0.07) x 10" 4 , (18) 

where the errors are statistical, experimental systematic, the estimated uncertainties from 
the niu form factor, and those from the piv form factors, respectively. 
For the decay B° — > p~l + v, we have determined the branching fractions 

B(0 < q 2 < 8 GeV 2 ) = (0.43 ± 0.20 ± 0.23 ± 0.09 ± 0.01) x 10" 4 (19) 
B(8 < q 2 < 16 GeV 2 ) = (1.24 ± 0.26 t° f 3 ± 0.22 ± 0.004) x 10~ 4 
B(q 2 > 16 GeV 2 ) = (0.50 ± 0.10 toil ± 0.19 ± 0.004) x 10~ 4 . 
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Combining these rates, again taking into account correlated systematic uncertainties, we 
obtain 

B(B° -> p~t + v) = (2.17 ± 0.34 ±° f 4 ± 0.41 ± 0.01) x 10" 4 , (20) 

where the errors are statistical, experimental systematic, the estimated uncertainties from 
the ptv form factors, and those from the ntv form factor, respectively. 

When the theoretical uncertainties that result from form-factor g 2 -dependence are evalu- 
ated in a common fashion, the branching fractions obtained in this analysis have uncertain- 
ties from the form-factor g 2 -dependence that are reduced by about a factor of two compared 
to previous ptv analyses [3, 4, 63]. These uncertainties are almost eliminated for the ntv 
branching fraction. 

We see evidence for the decay B + — > r]t + v with a statistical significance corresponding 
roughly to 3.2<r. The rate we obtain, 

B(B + -> rjt+v) = (0.84 ± 0.31 ± 0.16 ± 0.09) x 10~ 4 , (21) 

is consistent, within sizable errors, with that expected from the measured pion rate and 
isospin relations. Only an ISGW II form factor has been examined, and a 10% model 
dependence uncertainty has been assigned based on the previous CLEO ntv analysis. The 
final error quoted combines this estimate with the dependence on the irtv and ptv form 
factors. 

From the nlv q 2 behavior that we have observed, we find the ISGW II form factor for 
irtv consistent with data at only the 3% level. 

By fitting LQCD and LCSR calculations to the observed q 2 behavior in 7rtv, restricting 
each calculation to its valid q 2 range, and then combining the results, we extract 

\V ub \ = (3.24 ± 0.22 ± 0.13 ± 0.09) x 10~ 3 , (22) 

where the errors are statistical, experimental systematic, the estimated uncertainties from 
the ntv form-factor shape and normalization, and those from the ptv form factors' shapes, 
respectively. From a similar analysis of the ptv mode, we obtain 

\V uh \ = (3.00 ± 0.21 ±° H ± 0.28) x lO" 3 . (23) 

The errors are statistical, experimental systematic, theoretical systematic based on the 
LQCD and LCSR uncertainties, and ptv form-factor shape uncertainty, respectively. In 
general, the ptv form-factor calculations did not agree as well with the observed ptv data 
as did the ntv form-factor calculations with the ntv data. 

Combining these two modes for an overall result from this analysis, we obtain 

\V ub \ = (3.17 ± 0.17 ±g;i? t° H ± 0.03) x 10~ 3 . (24) 

Given the manner with which the theoretical uncertainties have been estimated, the quoted 
values should be interpreted as being closer in spirit to "one standard deviation" than to 
"the allowed range" . 

These results trade off the potential statistical gain over the previous CLEO analyses in 
favor of relaxation of theoretical constraints. Had we fixed the relative rate in the three q 2 
intervals in the ntv and ptv modes, a more pronounced improvement in statistical precision 
would have resulted. By relaxing the constraint, on the other hand, we have minimized our 
reliance on modeling in extraction of rates and of \V u b\. 
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These results supersede the iriv and plv results obtained in reference [3]. They agree, 
within measurement uncertainties, with the CLEO 2000 piv result [4] and with the recent 
BABAR piv analysis [63]. 

The results for \V U \\ obtained here are compatible with the results obtained from the 
recent CLEO end-point measurement [57]. The estimated theoretical uncertainties remain 
sizable for both itlv and piv, and there remain uncertainties in the estimates themselves. We 
therefore do not average these results, but view the compatibility as an indication that the 
uncertainties have not been appreciably underestimated. Significant progress in extraction 
of \V u b\ from exclusive decays will require a major improvement in theory. 

We thank A. Kronfeld, J. Simone, T. Onogi, T. Feldmann, P. Kroll, C. Maynard, and D. 
Melikhov for assistance with form factors. We gratefully acknowledge the effort of the CESR 
staff in providing us with excellent luminosity and running conditions. M. Selen thanks the 
Research Corporation, and A.H. Mahmood thanks the Texas Advanced Research Program. 
This work was supported by the National Science Foundation and the U.S. Department of 
Energy. 

APPENDIX A: DESCRIPTION OF EXPERIMENTAL SYSTEMATIC UNCER- 
TAINTY DETERMINATION 

The techniques employed in this analysis rest fundamentally on complete, accurate recon- 
struction of all particles from both B decays in an event. As a result, systematic uncertainty 
estimates that reflect uncertainties in the detector simulation must account for the reliability 
with which an entire event can be reconstructed, not just the signal particles. For example, 
if there is a residual uncertainty in the track reconstruction efficiency, the signal efficiency 
will not only be affected by incorrectly assessing the loss of the signal mode particles, it will 
also be affected by "misreconstruction" of the neutrino four-momentum. Furthermore, the 
rate at which background samples can smear into the signal region is also affected by the 
overall misreconstruction. 

We therefore estimate the systematic uncertainties due to detector modeling by modifying 
each reconstructed Monte Carlo event in each signal and background sample. For each study, 
the size of the variation has generally been determined by independent comparisons of data 
and Monte Carlo. The following list describes the variations that enter the systematic 
determination: 

tracking efficiency We have limited our uncertainty in track-finding efficiency for high 
(above 250 MeV/c) and low momentum tracks to be under 0.5% and 2.6%, respec- 
tively. These limits were obtained with hadronic samples, and therefore include any 
discrepancies in the interaction cross sections. To determine the systematic error from 
the uncertainty in tracking efficiency, we apply an additional inefficiency of 0.75% and 
2.6% to each high momentum track and to each low momentum track, respectively, in 
the simulation. 

tracking resolution We increase the mismeasurement of each momentum component for 
each reconstructed charged particle by 10% of itself, which is outside the range for 
which core distributions agree, but compensates for discrepancies in the tails. 

7 efficiency We have limited our uncertainty in photon reconstruction efficiency to 2%. In 
our studies, we have actually applied an additional 3% efficiency loss per photon, then 
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scaled the observed shifts back by 2/3. 

7 resolution We also degrade the photon energy resolution by 10% of itself. 

split-off simulation Studies of 77 — > K S K S have indicated that the combination of mis- 
modeling the physics processes and hadronic showers leads to an excess of isolated 
reconstructed showers (split offs) at the rate of 0.03/hadron in data relative to the 
Monte Carlo. To estimate the potential effect on our analysis, we interpret the entire 
excess as mismodeling of the hadronic showers, and add showers at this rate to each 
of our Monte Carlo samples. 

split-off rejection We bias our neural net parameter, which is derived from the distribution 
of energy within the crystals in the shower relative to the primary impact point of a 
"parent" charged hadron, to move photon-like results in the Monte Carlo towards 
hadronic-shower-like results. We limit the variations based on data and Monte Carlo 
comparisons of the parameter as a function of shower energy. 

K L showers In our simulation of K L showers, we increase the energy deposited in our Csl 
calorimeter. The variation is based on data and Monte Carlo comparisons of the energy 
deposited by showers after correction for the minimum-ionizing component. 

K L production By comparing the data and Monte Carlo K s energy spectrum and yield, 
we found that our Kl rate needed to be decreased by (7.2 ± 1.0)%, and that no 
correction was needed for the spectrum. The nominal analysis reweights events with 
Kl accordingly, and we vary the weight according to its uncertainty to estimate the 
systematic contribution. 

extra u production An important source of background is events that contain both a 
b — > civ decay and a c — > stv decay, where the latter can originate with either 
B meson in the event. We reweight the Monte Carlo so that the lepton momentum 
spectrum from secondary charm decay agrees with a spectrum obtained by convoluting 
a recent measurement of the charm meson momentum spectrum from B decay [64] 
with the MARK III measurement of the inclusive lepton momentum spectrum from 
charm decay [65]. The nominal result is corrected based on this procedure. To estimate 
the systematic uncertainty, we define spectrum "envelopes" and reweight our Monte 
Carlo samples to match this spectrum. The envelopes were defined by throwing 500 
toy Monte Carlo spectra in which all experimental inputs were varied according to their 
uncertainties and finding the variation within each momentum bin that contained 68% 
of the toy spectra. 

particle ID We simultaneously shift all dE/dx and time-of-flight distributions in the simu- 
lation by 1/4 and 1/2 of the intrinsic resolution, respectively. We take the full variation 
we observe as our uncertainty, even though this procedure leads to a very conservative 
systematic estimate. 

For each of these variations, we modify or reweight each event in each Monte Carlo sample 
in a full reanalysis of these samples. The set of modified samples for each variation replaces 
the nominal samples input to the branching fraction fit. For each variation, the shifts in the 
fit results provide the first input into the systematic estimates on the branching fractions for 
that variation. We can view the shifts in results as arising from two components: a change 
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TABLE XV: Central values and statistical uncertainties for B° — ► tt~1 + v branching fractions for 
the nominal fit and for each systematic variation of the Monte Carlo samples input to the fit. 
The detector-related systematic uncertainties in \V u b\ are obtained by fitting the results from the 
relevant set of q 2 intervals for each systematic study. The total branching fraction is shown as well 
for completeness. All results were obtained using the Ball'01 form factor for the irlv modes and 
the Ball'98 form factors for the ptv modes. 



systematic 10 4 x B(B° -► -k~1 + v) 

change total < q 2 < 8 GeV 2 8 < q 2 < 16 GeV 2 16 GeV 2 < q 2 < q 2 



nominal 


1 


.327 


± 


0. 


177 





.431 


± 


0. 


106 


0. 


.651 


± 


0, 


105 





.245 


± 


0. 


.094 


7 eff. 


1 


.348 


± 


0. 


194 





,476 


± 


0. 


117 


0. 


.674 


± 


0. 


117 


0. 


198 


± 


0. 


103 


7 resol. 


1 


.379 


± 


0. 


183 





.445 


± 


0. 


111 


0. 


.686 


± 


0. 


109 





.249 


± 


0. 


.096 


Kl shower 


1 


.311 


± 


0. 


173 





,426 


± 


0. 


104 


0, 


.642 


± 


0. 


104 


0, 


.242 


± 


0, 


091 


particle ID 


1 


.342 


± 


0. 


180 





,414 


± 


0. 


108 


0, 


.668 


± 


0. 


107 


0, 


.260 


± 


0. 


.096 


split-off rejection 


1 


.338 


± 


0. 


179 





,415 


± 


0. 


108 


0. 


.667 


± 


0. 


107 





.255 


± 


0. 


095 


track eff. 


1 


.357 


± 


0. 


185 





.446 


± 


0. 


112 





.669 


± 


0, 


110 


0. 


.242 


± 


0. 


.097 


track resol. 


1 


.317 


± 


0. 


179 





.438 


± 


0. 


108 


0, 


,664 


± 


0. 


108 


0. 


.215 


± 


0. 


.094 


split-off sim. 


1 


.326 


± 


0. 


178 





.432 


± 


0. 


108 


0, 


,655 


± 


0. 


106 


0. 


.240 


± 


0. 


093 


Kl production j 


1 


.325 


± 


0. 


176 





.431 


± 


0. 


106 


0, 


,651 


± 


0. 


105 


0. 


.244 


± 


0. 


.094 


Kl production { 


1 


.330 


± 


0. 


177 





.432 


± 


0. 


107 


0, 


,653 


± 


0, 


105 


0. 


246 


± 


0. 


.094 


v production f 


1 


.344 


± 


0. 


178 





.425 


± 


0. 


106 





,669 


± 


0. 


106 


0. 


251 


± 


0. 


.095 


v production \ 


1 


.322 


± 


0. 


175 





.439 


± 


0. 


106 


0, 


.641 


± 


0, 


104 





.242 


± 


0. 


.093 



in the signal efficiency and a change in the predicted background level. These changes tend 
to cancel in the total shift: a variation that reduces the signal reconstruction efficiency also 
simultaneously increases the background level (and reduces the signal yield from the fit). 
As the main text describes, we increase our systematic estimate to allow for imperfections 
in the predicted cancellation. 

APPENDIX B: EXTRACTION OF \V ub \ FROM THE MEASURED dT(B° -> 
tt~1+v) /dq 2 DATA WITH FUTURE FORM-FACTOR CALCULATIONS 

The branching fractions in the three q 2 ranges for B — > niu exhibit very little dependence 
on the precise form factors used to extract the branching fractions. The results can therefore 
be reliably used to obtain values for \V u b\ using future B — > niu form- factor calculations 
that are improved over those used in this paper. This appendix provides the detail needed 
to ascertain the proper experimental uncertainties for such an extraction using the same 
fitting technique presented above. The main difficulty stems from proper evaluation of 
the experimental uncertainties because of correlations (both positive and negative) among 
the results for the three ranges. The correlations arise both statistically from the fitting 
procedure used to extract the three rates and systematically as we vary the details of the 
simulation. 

To extract a central value of \V u b\, we perform a \ 2 fit to the nominal branching frac- 
tions from the three q 2 intervals listed in Table XV. This \V u b\ fit includes the correlation 
coefficients among the rates from the branching fraction fit to the data: p i2 = —0.035, 
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TABLE XVI: Fractional uncertainties to be added in quadrature to systematic shifts in \V u b\ to 
account for uncertainty in cancellations arising from correlated efficiency and background changes. 
The correction is shown for the various different q 2 ranges used in this analysis. 



systematic 






additional systematic i 


;%) 




change 


full range 


< q 2 < 16GeV 2 


< q 2 < 8 GeV 2 8 < q 2 


< 16 GeV 2 


16 GeV 2 < q- 


7 eff. 


1.67 


0.51 


0.72 


1.22 


1.49 


7 resol. 


0.19 


0.28 


0.14 


0.43 


0.30 


Kl shower 


0.25 


0.30 


0.14 


0.16 


0.46 


particle ID 


0.25 


1.09 


0.29 


0.27 


0.58 


split-off rejection 


0.00 


0.56 


0.24 


0.21 


0.35 


track eff. 


0.99 


1.62 


0.72 


0.90 


1.17 


track resol. 


0.49 


0.25 


0.14 


0.11 


0.44 


split-off sim. 


0.23 


0.39 


0.24 


0.11 


0.17 


Kl production 


0.01 


0.02 


0.01 


0.00 


0.01 


v production 


0.12 


0.43 


0.28 


0.19 


0.13 



p l3 = 0.003, and p 23 = -0.037. 

To evaluate the error arising from simulation uncertainties ( 'V simulation" in Table IV) 
on the results, we redo our x 2 fit for \V u b\ using the new rates listed in Table XV for each 
variation. For the results presented here, we have used the correlation coefficients from the 
branching fraction fit to the data for each variation. In practice, the coefficients remain 
stable enough that using the nominal coefficients in all fits is sufficient. The change relative 
to the nominal \V u b\ result provides the first input to the uncertainty estimate. For the 
uncertainty estimate in K L production and secondary v production, we take the average of 
the "up" and "down" shifts as our overall estimate. To allow for misestimation of correlated 
changes between background levels and signal efficiencies in the results (see main text), we 
increase the fractional uncertainty on \y u \\ from each variation by adding in quadrature the 
quantities listed in Table XVI. Finally, the 7 efficiency uncertainty should be scaled back 
to 2/3 of the value found above. We combine all of the uncertainties in quadrature to arrive 
at the total 'V simulation" systematic for \V u b\. 

We evaluate the uncertainty from our modeling of the B — > X u tv backgrounds in much the 
same fashion. The fit variations that we have used for this purpose are listed in Table XVII. 
An earlier version of our B — > X u tv generator was used in the study, and the table also 
shows the "nominal" result obtained with that version. We did not expect large differences 
from our change, and indeed the results obtained are very similar to the nominal results 
in Table XV. To obtain the uncertainty estimate resulting from the hadronization model, 
we compare the results using purely nonresonant hadronization to that using our nominal 
mixture of resonant and nonresonant modes. To obtain the uncertainty resulting from our 
choice of parameters for the OPE-based inclusive differential rate calculation, we take the 
average of the shift from the last two lines in the table relative to the above nonresonant 
result. Note that these variations do not affect our signal Monte Carlo samples. 

For the remainder of the systematic uncertainties, we take one half of the fractional 
uncertainties listed in Table IV. The factor of one half arises because of the square root 
involved in extraction of \V u b\ from the rates. 
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TABLE XVII: Central values and statistical uncertainties for B° — > tt~1 + v branching fractions for 
the reference fit and for each systematic variation of the B — > X u lv background simulation input 
to the fit. The associated systematic uncertainties in \ V U \\ are obtained by fitting the results from 
the relevant set of q 2 intervals for each systematic study. The total branching fraction is shown as 
well for completeness. All results were obtained using the Ball'01 form factor for the ttIv modes 
and the Ball'98 form factors for the ptv modes. 



OPE 


hadron- 




10 4 x 


B(B° -► n~t+v) 




parameters 


ization 


total 


< q 2 < 8 GeV : 


2 8 < q 2 < 16GeV 2 


16GeV 2 <g 2 <g max 


nominal 


nominal 1.324 ± 0.177 


0.423 ±0.107 


0.655 ±0.105 


0.246 ± 0.094 


nominal 


nonres. 


1.322 ±0.177 


0.431 ±0.106 


0.639 ±0.105 


0.251 ±0.094 


"High" 


nonres. 


1.311 ±0.176 


0.428 ±0.106 


0.637 ±0.105 


0.246 ± 0.094 


"Low" 


nonres. 


1.329 ±0.177 


0.434 ±0.106 


0.646 ±0.105 


0.248 ± 0.095 



APPENDIX C: AVERAGING \V ub \ RESULTS 

In each of the ntv and ptv modes, we have extracted two results for \V u b\ that are largely 
free from modeling assumptions: a value based on the application of LCSR-derived form 
factors for q 2 < 16 GeV 2 , and a value based on the application of LQCD-derived form factors 
for q 2 > 16 GeV 2 . We therefore have three averages to be calculated: the combination of the 
two results within the ntv mode and within the ptv mode, and the combination of the two 
modes. The averaging procedure should take into account, in particular, the correlations 
present in the systematic uncertainties in the result. This appendix describes our averaging 
procedure. 

The statistical correlations have been taken into account in the LCSR-derived results. An 
evaluation of remaining statistical correlations found that they had little impact on the final 
statistical error, and we have not included them in the final procedure. Proper treatment 
would have led to a decrease in the overall uncertainty that would be hidden at the quoted 
precision. 

Regarding theoretical uncertainties, while the two techniques have different systematic 
effects, both approaches currently have systematic issues that are difficult to evaluate. For 
example, there is a quark-hadron duality assumption in the LCSR approach, and the current 
LQCD calculations have been evaluated in the "quenched" approximation. Treating the un- 
certainties as uncorrelated would therefore be likely to underestimate the "true" theoretical 
uncertainty. To be conservative, we treat the theoretical uncertainties as if they were fully 
correlated. 

Let us first consider the two results obtained within a given mode. We wish to combine 
the results with a weight that minimizes the overall uncertainty and preserves the systematic 
correlation information. Defining the weight of the LQCD-derived result (denoted ll^f,]- 16 ) 
by a, the LCSR-derived result (denoted |V^6| <16 ) enters with a weight 1 — a: 

\V ub \ a = a\V ub f 16 + (1 - «)|K6| <16 . (CI) 
The statistical uncertainties are uncorrelated, and are combined as 

4at = («^at 6 ) 2 + ((1 " «R1at 6 ) 2 - (C2) 

Correlated uncertainties, such as the theoretical uncertainties, are combined as 

<W = («^ rr 6 ) + ((1 - a)a^ T ). (C3) 
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For each simulation variation (labelled i), we perform the full analysis to obtain |14&L^ 16 ) 
and \V ub \f 16 . The systematic uncertainty defined for the variation is 

<n = \V ub \7 m - («|K4- 16 + (1 - u)\V ub \< 16 ), (C4) 

where |14&la° m is the average resulting from Equation CI. This procedure preserves the 
systematic correlation. We combine this estimate in quadrature with the additional uncer- 
tainty contribution to allow for imperfect modeling of the correlated changes between signal 
efficiency and raw yield (see Section V). 

Finally, for each value of a the experimental and theoretical uncertainties are combined 
in quadrature (taking the average theoretical uncertainty in the case of asymmetric uncer- 
tainties). We scan over a and choose the value that minimizes the total uncertainty. 

We perform a similar procedure to combine the results from the two modes. The weights 
obtained individually for the different q 2 regions in each mode are fixed. The uncorrelated, 
correlated, and anticorrelated uncertainties are combined in exact analogy to the above 
descriptions. Taking (3 as the weight of the iriv mode in the average, we have 

\V ub \p = f3\V ub Y + {l-f3)\V ub \ p . (C5) 

For each simulation variation, the systematic estimate becomes 

<Ji = \V ub \f m -[p(a n \V ub \f 1 ^ + (1 - a n )\V ub \f w n + (1 - (3)(a p \V ub \f w > p + (1 - a p )\V ub \f 16 '") 

(C6) 

These uncertainties are, as before, combined in quadrature, along with the contribution for 
imperfect modeling of the correlated efficiency and yield changes. 

We scan over the weight (3 to find the value that minimizes the overall combined uncer- 
tainty. Once again we treat the theoretical uncertainties in the ntv and plv form factors as 
correlated in this procedure. 
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